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Sommaire. — La thermodynamique des processus irréversibles, qui repose sur l’emploi des relations 
réciproques d’Onsager, est développée pour un mélange dans le cas de la diffusion thermique (effet Soret) 
et de phénomènes connexes comme l’effet Dufour (effet thermique de la diflusion) et les potentiels de 


Pa diffusion ordinaire et thermique. 


# 


1. Introduction. — La thermodynamique des 


processus irréversibles repose sur l'emploi des 
_ relations réciproques d’Onsager [1]. Dans ce mémoire 
_ une théorie sera développée pour un mélange de 
corps chimiques dans lequel s’opèrent les phéno- 
_mènes irréversibles de conduction de chaleur, de 
_ diffusion et de conduction électrique ainsi que les 
effets résultant de l'interaction de deux d’entre 
eux. La diffusion thermique par exemple (qu’on 
appelle effet Soret [2] quand il se produit dans 
l'état condensé), se manifestant par la diffusion 
de matière dans un champ de température non 
uniforme, résulte de l'interaction de la conduction 
- de la chaleur et de la diffusion. Il existe pourtant 
. un autre effet dépendant de cette même interaction 
_et qui est en un certain sens l’effet inverse de l'eftet 
Soret, c’est l’effet Dufour [3] selon lequel un gradient 
de température s'établit quand deux corps diffusent 
l’un dans l’autre. Or, d’après la théorie d’Onsager, 
-il y a des relations réciproques entre phénomènes 
irréversibles inverses. l > 

Pour les phénomènes traités dans ce mémoir 
_ il existe aussi des théories « quasi thermostatiques », 
dans lesquelles on applique les lois de la thermo- 
- statique habituelle à une partie du processus sup- 
‘posée réversible. Eastman [4] et Wagner [5] ont 
_ développé une pareille méthode pour l’état station- 
naire de la diffusion thermique. Ces auteurs émettent 
hypothèse qu'un transport virtuel de matière 


peut être considéré comme réversible et ils font 
abstraction de la conduction de chaleur irréver- 
sible qui subsiste à l’état stationnaire : la méthode 
est évidemment erronée, mais les résultats sont 
confirmés par l’expérience et par la théorie thermo- 
dynamique (c’est-à-dire avec emploi des relations 
d'Onsager). On peut également construire une 
théorie « quasi thermostatique » pour les potentiels 
électriques dans le cas général non stationnaire [2]. 

Dans ce mémoire la théorie thermodynamique 
avec relations d’Onsager est développée dans le 
cas de la diffusion thermique et les phénomènes 


_ connexes, et l’on trouve entre autres résultats toutes 


les formules de Wagner pour l’état stationnaire 
et les expressions pour les potentiels électriques. 


2. Les paramètres conjugués. — Pour trouver 
les paramètres conjugués, qui entrent comme 


« flux » et « forces » dans les relations phénomé- 
nologiques, décrivant les processus irréversibles, 


il faut calculer la source d’entropie, quand il y 


a des fluctuations dans le système considéré. Je 
veux généraliser le raisonnement d’un premier 
mémoire [1] dans deux directions. Premièrement 
j'admettrai que le système est un mélange de } 
composants (1 — 1, ...,j), et deuxièmement je tien- 
drai compte de l'existence éventuelle de forces 
externes Æ; sur l’unité de masse du composant à 
afin de pouvoir déduire les effets électriques. 
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Supposons que le mélange se trouve dans un 
récipient à deux parties. Alors on peut écrire pour 
la variation de l’entropie dans la première partie, 


à cause des fluctuations AU dans l’énergie et A, m 
dans la masse du composant à 


S 
== = PAU+ > 9 D Ain + LU 
2ŒS 
ce à on AU A mi 


Fe Don Fa _. AimiAime (4) 


et pour la seconde partie une expression semblable 
avec indices 2. Supposons que le récipient soit 
adiabatiquement isolé. La loi de conservation de 
masse s'exprime comme 


Aimi+ Aom;= 0 (A RESTE nie (2) 
et celle de l’énergie comme 
NU AA DEN Am; Ax, (3) 


où Ax est la distance sur laquelle la masse Am; 


est transportée. Si l’on écrit 
2AS = A S+AS (4) 


pour la variation de l’entropie de la masse 


2 D mi = 9 M; 


on trouve 


AS = FN Am Ar + = PAU) 


ds 


JU dm SU Ami + = Do ve (5) 
à, k 


si l'on n’emploie de ee que jusqu’au ue 
ordre et avec 


AURDAUS—AU; et Ami= im 


La dérivée temporelle de AS est égale à 
k pe ë ES 222$ 
HAE DE 


i 
22, 
TAN Cor se 
S _2S Ne 
F2 om; mx UGS Te Fi ré 2 


dS 
NS AU 1a(S) 


. JS DÈnEe | Fa 
+Dam af) D me Jo lit f Qu) 
i 
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Si (7) est dort sous. dore ce somme de rod 
de paramètres conjugués 


AT (4 e +Y Te ) Ar, 


les « flux » suivants sont introduits 


Jy=AU (courant d'énergie) (95% 
it 2 
J;=Am ( » de matière é). (10). 
On a les relations thermostatiques CR. 
NS I 2 
—— —= — 11% 
oU 1 ss ) 
et SSSR 
OS. pu 19 
Ë ANA ENT ce 


de sorte que les « forces » correspondant aux « flux » - 
sont. 


: AT AT grad 7° LA 
UP aTT ENT 2 ce 
et 
A NE Aie ’ 
A a(#) re sad (5) + 7? (14) 


avec p le potentiel chimique du composant 1. 

La « force » (14) est composée de deux termes. 
Le premier [6] n’est pas simplement le gradient 
de sn comme dans la loi de Fick, mais 
le gradient de - Le second terme contient la force 
externe, ici la Fe électrique 


(5) 


où e; est la charge de l'unité de masse du composant à 
et où © est le potentiel électrique. Le courant élec- 
trique J; est relié au transport de matière de la 


F;=—e;grado, 


substance à i par 


Ti ="eil;, 
de sorte que courant électrique total est 


(l SJ ei 
i i 


La contribution du courant électrique à AS contient 


VE 
> JiFi== > T; grado = — J grade, 
Où i 


= 


(18) 


la chaleur de Joule. 
Comme le potentiel chimique est une fonction 


de la température, de la pression et de la compo- 


sition, exprimée par (j—-1) des j concentrations: 


relatives n; ps que > V i=1 an] à 


L 


= n L à ù . 
Di (I, p, nn, .., Riu), 


es 


sh panda in at aan dpi ta 


V 


«Srad ut; — = surad7 +oigap eY 2 = grad nx 


= 


; 3 sur 
LE VAT Cal Pa VEUT re 


RE age 


(21) 


KT 


{s: et v;, entropie et volume molaire partiel du 
composant < se sorte que la force X; peut s’écrire 


16e TE grad 7° »; gradp 


Fe 
2 


ni 


du; 
(SE : :) grad nx— €; grid, (22) 


où w, est Périhalpie molaire partielle du composant 


( CRETE ES Dsy): 


3. Les relations phénoménologiques. -- Pour 
un mélange isotrope on peut admettre comme 
première approximation 


à Ji =Y ax Vi+ DiXu 1 1 Er CD 
. k 
Ju=Ù br Ar + cAu. (2%) 


k 


On voit dans (23) que les forces X4 donnent 


naissance entre autres choses à des flux aXx. 


Or, 


suivant (14) ou (22), la force contient des 


termes en grad n4 et grad +. Dire que ces « forces » 
engendrent des flux qui leur sont proportionnels 


c'est exprimer les lois de la diffusion et d’Ohm. 
Pour les systèmes isothermes et idéaux la loi de 
la diffusion est identique à la loi de Fick, qui dit 
que le flux de la matière est proportionnel au gradient 
de concentration. Il résulte de (24) qu'un gradient 


de température T°? X, provoque un courant d’énergie 


cX,, expression de la loi de Fourier. On voit éga- 
lernent se présenter dans les relations phénomé- 


. -nologiques des termes provenant d’une «interaction » 
de deux phénomènes. Par exemple azX% est le 


flux de matière de l’espèce 1 résultant d’une force X4, 
liée à la substance k. Le dernier terme à droite 
de (23) est un flux de matière b;X, provenant d’un 
radient de température T?X,. Il représente donc 


une diffusion thermique. Ce n’est toutefois pas la 
diffusion thermique entière parce que X; contient 


lui aussi un terme avec grad 7. Dans (24) on voit 


que le terme b;X; contient l'effet inverse de la 


- diffusion thermique 


c’est l’effet Dufour, celui 
d’un flux d’énergie provoqué par un gradient de 


_ concentration. 


4. Les relations d’Onsager. — La première 


classe de relations réciproques d'Onsager s'exprime 


par 
ik = ki : (25) 


et montre que le flux de matière de l'espèce 1 causé 
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par une « force » X,4 est régi par un coefficient qui 
est égal à celui qui gouverne le flux de k causé 
par X;. 

Plus intéressante encore est Ja deuxième classe 


de relations d’Onsager 
Dy= D}, (26) 


reliant les coefficients de diffusion thermique et 


‘ d'effet Dufour. 


D. L'état stationnaire. -- Introduisons d’abord 
les grandeurs Q; définies par 


= aix Q4 CRT D (27) 
Æ 
On obtient alors en introduisant (27) dans (23) 
= auri+ OPEL USA, 2 UNIT C8) 
' 
La définition de l’état stationnaire est 
OMG RO mais J, < 0: (29) 


Les relations (28) deviennent donc pour ce cas 


A+ Qi Au 0 CRC 5 (30) 
ou, en écrivant explicitement les forces à l’aide 
de (13), (14) et (15) 


grad u;+ 6e; gradg +(OQ;— 1) T1 gradT = 0 


o 
Vie NS, ND) 
où encore avec l’expression (21) ou (22) 
fal 
du; 
e; grad p +Y ne grad n# 
RE: À (32) 
+ e; grad 9 + (0; w;) 171 grad7T'= 0 


(i=T, VD 


Telles sont les équations générales régissant l’état 
stationnaire. 

Dans la théorie présente les relations d’Onsager 
n'ont pas encore été utilisées. On a vu d’autre 
part que les grandeurs Q; n’ont été définies que 
formellement par (27) comme des combinaisons 
des coefficients a; et b;. Or, on peut démontrer 
qu'il est possible de donner une interprétation 
physique des Q; à l’aide des relations d’Onsager (25) 
et (26). À cette fin multiplions les équations (23) 


par Q; et additionnons-les 


NT — à Fe 
D Qi= D Qian let O6: LE (33) 
l ih L 


En soustrayant cette équation de (24) on trouve 
que les premiers termes des seconds membres 
s’annulent comme conséquence des relations d’On- 
(25) et. (26) appliqués. & (27). On a par 


- s ER ne ue 
; : € ss RE % $ ; A 
LAS re ES ET 
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conséquent En multipliant les ne (32) respectivemen 


(34) 


7 =D + or ÿr. 


Cette équation permet de donner une interpré- 
tation physique : à l’état isotherme (X, — 0) 
Q; est l’énergie transportée par l’unité de masse 
du composant à (c’est-à-dire quand J;— 1). On 
pourrait appeler Q; une « énergie ou chaleur de 
transport », si cette expression n'était déjà utilisée 
pour la notion 


= Q;— 3, (35) 


où Q; est la chaleur absorbée par une partie I 
du système (avec température 7) et dégagée par 
la partie II (avec température T + dT) quand 
. l'unité de masse du composant i passe de I à IT, 
de sorte que les températures de I et II conservent 
leurs valeurs initiales, la pression étant uniforme. 

La formule 60 montre que pour l’état station 


naire ; 
10e a > Qib D 7e) (36) 


tandis qu'aussitôt après l'application du champ 
de température non uniforme, c’est-à-dire quand 
“on a encore grad ny = o (k = 1,'..., j), geo 


et grado = 0 
“Ju Ë re) je 


pi 


(37) 


comme on voit en introduisant (22) dans (24). 

Le coefficient de conductibilité de (36) contient 
les grandeurs Q;, qui sont caractéristiques pour la 
migration de matière qui a eu lieu pour arriver à 
l’état stationnaire. Le coeflicient de (37) contient 
au lieu des Q,; les grandeurs thermostatiques w:. 


6. L'état stationnaire. La diffusion thermique. 
— On dispose maintenant comme équations géné- 
rales de l’état stationnaire de (31) ou (32) avec, 
si l’on préfère employer Q; au lieu LE Q:, les rela- 
tions (35). 

Les équations (32) avec 


CT Q | 
> ER ou > gradn;= 0 (58) 


i 1 U 


et la condition de la neutralité électrique 
D'uer ou Da 


forment un système de j + 2 équations linéaires 
qui permettent de calculer les grad ny (k = 1, …., j), 
- grad® et grad p comme fonctions de grad T. Je 
veux ici déduire les équations différentielles de 
la diffusion thermique (ou effet Soret, quand il 
se produit dans un liquide ou un solide). 


grad;= 0 (39) 


Par Ni, Nos cs Nj et en les additionnant, on trouve, 
eu égard à la relation de Duhem-Margules dr dui= ra 
et comme conséquence de (39) 


PA grad 7 
p grad p + (Ze) D = 0 


i 


pr > AjVyr. 
è 


Dans les problèmes à envisager ici on 
toujours que la pression est tonne: 
de (40) qu’alors à 


d'aiQ 0: 


i 


co) # 
avec 


observe … 
on tire 


| a 


Les chaleurs ae transport Q ne peuvent donc pas Ë 
avoir toutes le même signe. É 

Trouver les équations différentielles de la diffusion 3 
thermique revient à trouver les grad n; des équa- 
tions générales (32), (38) et (39). Ici, quelques 
cas particuliers seront considérés sur les mélanges # 
isobares et idéaux. Les qu (32) s écrivent | 
dans ce cas { 


crad »; : opad hs ET 
gr rad * + e;grado + O2 = (4) 


R F6) 
ni À s T : 
parce que 
Ou; LS RT 
ARR RES 


: (43). | 


où les substances à — x, , (j—1) sont les «corps 4 
dissous », tandis que j est le « solvant ». On trouve 
ce que l'on appelle le « coefficient de Soret » si l’on. 
élimine grade. A cet effet il faut multiplier les 
équations (42) respectivement par me; et addi- 
tionner É à e S 


ni QŸ $ S 
i grad FAi (44) À 


e 3 ÿ ° : 
CPI = a 


i : ! LA 


En introduisant grade de (44) dans (42), on Ë 
trouve pour le coefficient de Soret L 


_ 
€; der nr Of 


I grad”; Of Fr 
ni SrAd EE: REF LT 
Fe, RTI etnx 


k 


Pour un composant non chargé, €; = 0, le dernier. 
terme est nul 
I Era à SENOT 


ni gradT 


Si l’on s'intéresse à la séparation relative de deux. 
composants à et X sans charges électriques, deux. 


fes emple, il faut connaître Ta ture 


e : n; 


grad j 
RE . grad À 
I} 
Fr grad”; e 1. gradnx _ Qj— Of ES 
= An; grad 7° ñn}x grad 7° + DEPTe (4 


Pour la solution d’un électrolyte avec cations 
et anions monovalents (e—:, @—-6, € charge 
élémentaire de l'électricité) dans un solvant neutre 
(e — o) l'équation (45) devient 


1 gradz: HN grade 
no grad 7' 


ni atad TS 0 
2 
= GA ÆQ HRQE" 
Rte es ARTE 3 RT? ce 


Un cas intéressant qui ne peut se réaliser que dans 
une solution solide est celui d’un sel contenant 
deux cations « 1 » et « 2 » mobiles et un réseau 
fixe d’anions « 3 » (grad n; — o). Alors on trouve 
pour l'effet Soret des cations dans une solution 
_ idéale 
grad 7;  gradñns 
7 gradT grad7' 
Qi—Qi _ Qi—Oi rm 
r(5 +) Le RT? MEne 
on: on 


2 


(49) 


Avec une notation plus commode, employant 
_ les fractions moléculaires des sels 1 3 et 2 —3 : 


7; UE) 
EU re 5, à () 
Hi + is Ii + 79 


il résulte que 


- 5 : 1 gradni3 I grad no; 
| FE grad 7 F es grad 7" 
I grad 722 Qi — Q: 
Re moderne po 
On a vu que l'effet Soret dépend des chaleurs 
de transport Qf = Q;—w;. Il est nul quand ces 
chaleurs disparaissent ou quand Q;—w;. On voit 
aussi que dans ce cas les coefficients de conduc- 
 tibilité de la chaleur de (36) et (37) sont égaux. 
Cela est évident puisque dans un mélange qui ne 
présente pas d'effet Soret, il ne s'opère pas de 
migrations de matière et la conductibilité de Ja 
chaleur est donc toujours la même. 


7. L'état stationnaire. L'effet thermoélectrique 
_ homogène. — Les équations générales données 
au commencement du paragraphe précédent montrent 
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qu’ un gradient de potentiel s’établit dans le mélange. 
‘Ge gradient de potentiel, qui est trouvé en résolvant 


grad 9, est appelé l'effet thermoélectrique homogène 
parce qu'il se trouve dans l’enceinte du mélange, 

Considérons d’abord les mêmes exemples que 
ceux traités dans le paragraphe précédent, Pour 
un mélange idéal on a trouvé l’équation (44). 

Pour une solution d’un électrolyte avec e, — € 
et &——e dans un milieu électriquement neutre, 
il s'ensuit donc 


Il u 
grad o =-— = (Qi — Q* ee (32) 


Dans le dernier exemple traité dans le para- 
graphe 6, on a grad © — 0. 

Il existe toutefois un cas spécial qui est réalisé dans 
certains électrolytes solides. C’est celui où l’un des 
deux composants est mobile et où l’autre forme le 
réseau fixe (exemples : électrons dans un métal, 
cations mobiles d’une structure « 4 » avec anions 
fixes). Alors l'effet Soret est pratiquement impossible 
à cause des grandes forces de Coulomb : grad n;, = 0. 
L'effet thermoélectrique homogène est donc (42) 


Of grad7 
re i Sladz Wa 
Érp er  Emt (53) 
2j 1 


8. L'état non stationnaire. -— Dans l’étude 
de l’état stationnaire (J; — o) un potentiel élec- 
trique qui se présente dans le mélange anisotherme 
a été trouvé ($ 7). On peut se demander quelle 
est la valeur du potentiel électrique quand l’état 
stationnaire n’est pas encore atteint, en d’autres 
mots quand on a pour les flux de matières J, 0, 
(—1, ...,j). À cause de la grande interaction 
coulombienne le courant électrique total est toujours 
pratiquement nul 


2e ve ETS (54) 


i 


c’est-à-dire que le mélange reste pratiquement neutre 
électriquement. 


En combinant l'équation (54) avec l'expression 
de Je (28) 


li =Y ax Ar+ Qx A), (55) 
ke 
on déduit 


Dead ant Fe+ Qi du) = 0: (56) 
er ke 


a écrivant explicitement X, (13) et X; (14, 15) 


on obtient 


pee je 
pics) Gite (- F grad =. 
è k 


grad 7° LEA 
— ex grad 9 — Or —— ins )= 0: 2001) 


NE ” 
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pres #1 ; 2 de de 
donc TN 
a! | Lx = 
Daire; grad À d'ane:Qu 
me mn HAUTES 
grad @. — = — - = - ABRS (58) 
TN aire;e, djpe;er 
dd 
ik i,k Ç 
_ Avec l’abréviation 
ze 
> aikeiek 
RESUME A (59) 
> direiek 
: ik 
on en tire 
UE L' l orad 7 
grade = T'erad(#à) >. Où > 7 (60) 
k Fe : : f 
ke be grad 7 
grade =—Ù * erad u4 > SHOP CREER GI 
5 è er? PA ex CC be) T K ) 
k s Æ 
Les relations d'Onsager (ax— ax) permettent 


d'interpréter les grandeurs f; (59) comme « nombres 
de transport » {; du composant k. On entend par 
nombre de transport {; la contribution que le 
corps kX fournit au transfert d'électricité. Afin 
de démontrer ceci, commençons par écrire le courant 
électrique que le composant k produit . quand 
grad T = 0 et quand la « force » X;—=e; T-!grado 
n’est causée que par le champ externe grad () 


1x = er = er N a; V;—=— 0e} ape; T1 grad o. 
LR kJK r> fi Ai RS krCi S 


u i 


(62) 


Le nombre de transport a est défini plus haut 
comme le quotient 
> akie;ek l 
Yr Æ à dkiC;ek 
k 


à 


N 


Si l’on applique la première classé des relations 
d’Onsager ax = &x, on trouve que 


le = Uk, 


(64) 


de sorte que les accents dans les équations peuvent 
être supprimés. 


9. L'état non stationnaire. L'effet thermo- 
électrique homogène : potentiel de diffusion 
et potentiel de diffusion thermique. a- 
tion générale du potentiel (61) peut être récrite 
si l’on y introduit l'expression de gradp; (21) 
et Q; (35). 


On trouve 


! 


grad 7 
T 


os sy (88) 


grade Ù À Of 
€ 


_ 


de son contenu physique. 


‘aux coefficients de diffusion ordinaire D. En effet, 


RE E VANCET 
ee = L Ê 
ci # à es at ( 
S LIST 2 EE 


< En s a 


celui dé mélange deal 1743) dt ban HR 


“ire 


A 3 
S Aiie( 
i 


Ar = Ôjk ji, donc 7; = 


Alors l'équation générale (65) est: 


gradoy—+ grad 2 3 4 
di a Dir (67) | 


€i 
i 


grad © — 


Cette équation se prête très bien à une discussion | 4 


On peut naturellement trouver directement Ja 
relation (67) en appliquant (54) à l’équation du. 
flux de matière (13, 22, 28, 35) dans un mélange 
idéal et isobare (66) AE 


Gi gra d 7; ee ae 
Je a (- RIT à — e; grade — Q; }. 


(68) 4 
= 

Le premier terme du second membre de (67), 
grades, est l'expression connue du potentiel de . 
diffusion ordinaire qui apparaît quand les vitesses 
de diffusion des ions sont différentes. Ce terme. 
contient les nombres de transport {; qui sont reliés. 


D; est par définition le coefficient de gradn, do À 
l’expression du flux J; (68). Par suite 


&;; R 


D;= L Qer 
et avec (66) 
a n; D;e : (70). 
dniDie : 
: = 
J= Di(— grad 7; — SEP grade — QUE rad r. (70) 
RATS ns | 


La relation (69) est identique au théorème de 


Fokker-Einstein qu’on écrit souvent ‘ 

» Lt à 

D;= U;RT, 118 

avec U, la « mobilité » de à. Cette mobilité est la ” 

vitesse moyenne de À quand la force extérne X;= 1 Ë 

(73) | 

Pour X;=1 on a nU;=— «ax et donc (69), comme il 
convient. 

L'autre terme de (67), grado,, traduit l'influence 
de la non-uniformité de la température. Le cosffi- … 
cient de gradT' contient les produits des nombres 
de transport /; et des chaleurs de transport Q. 
Ces produits sont une mesure de la diffusion ther- « 
mique. Cette deuxième contribution à grado existe | 
donc quand les chaleurs de transport des ions sont 
différentes. Il paraît donc. naturel d'appeler ce. À 
terme «le potentiel de USION thermique » PR 


Ji=nioi= Qi A. 


inalogie ave 


y 


. trouvé pour l'état non stationnaire, c’est-à-dire 


à un instant arbitraire. Considérons maintenant 


_ comment grade varie dans un mélange originai- 


rement homogène qui est plongé dans un champ 
de température gradT au moment {—0o et dans 
lequel s’opèrent les phénomènes de diffusion ther- 
mique et de diffusion ordinaire jusqu'à ce que 


_ l’état stationnaire soit atteint ({— æ). Au moment 


où {— 0 les concentrations sont encore uniformes 
(gradn;— o), de sorte qu’il n'existe pas de potentiel 
de diffusion ordinaire. On a donc 


N' 40 grad7 
Lil ei T° 


è 


grad 96 = — 


; (74) 


c'est-à-dire que grade pour {—o n'est que « le 
potentiel de diffusion thermique ». Pour la solution 
d’un sel avec e, — s et e, = (74) est 


grad 7° 


Seralqu = (A Qi Qi) 


(75) 
Ce potentiel peut être déterminé aussitôt après 
le commencement d’une expérience. 


Au fur et à mesure que les gradients de concen- 


. tration s’établissent à cause de la diffusion thermique, 


un potentiel de diffusion se fera sentir. Le potentiel 


- total est donné par (67). Pour la solution d’un 


électrolyte avec e, — € et e —— € 
grado = gradsy+ grad 6, (76)' 
RT. grad 7 
grade = (at) 
, srad 7° Le 
(4 QT — f05)° : (77) 


A l'état stationnaire (f— æ) on a Ji—o et 
pour l'effet Soret on trouve (48). En introduisant 
cette valeur dans l'équation générale de grado 
on retrouve pour le mélange idéal 


DT 


grad 7” 


i 5 
PORTE ARE TRS 18 
sradé, Ee T- (78) 
ve 1; : 
s L 
et pour e, — € et ee ——E€ 
do Di OS grad7 
DC DT d'os ei Ca Ti 
REZ -. grad 7 
—-(0,.--0;)> 16 
60: 23) T u : (79) 


_ équations déjà déduites [(44) et (52), et qui ne 


contiennent plus les nombres de transport 4. Cela 
est évident puisque à l’état stationnaire le transport 
est nul. 


| TR à RRUpE Cie # ï à ee : h ae 2. L à me 
El MOI YNAMK UE DE QUELQUES PROCESSUS IRRÉVERSIBLES 
alogie avec le premier terme qui est l'expression 
du potentiel de diffusion ordinaire. ie : 
L'effet _thermoélectrique homogène (67) a été 


= 


/ 
f 
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Dans la figure 1 l’effet thermoélectrique homogène 
dans une solution aqueuse d’un électrolyte avec 
ee et &=—-e est porté en fonction du temps. 


8| 107 Voit/cm 


_—— 
en 


{ à 
20 heures 


M T6 1 


Fig. 1. — If*effet thermoélectrique homogène (grad+) comme 
somme du potentiel de diffusion thermique (gradvw,) et de 
diffusion ordinaire (gradvw,) en fonction du temps b). 


Pour calculer cette fonction il faut savoir comment 
grad n varie. Cela ressort de la théorie phénomé- 
nologique [7] qui dit qu’approximativement grad n : 
s'établit suivant une loi exponentielle 


L 
grad nr “ 0 5) 


(80) 


a? à k Re 
pp Où a est la hauteur d’un récipient 
dont la paroi supérieure a une température plus 
haute que la paroi inférieure. D est le coefficient 


de diffusion. Pour le cas envisagé on peut écrire (80) 


avec OÔ — 


- en vertu de (48), qui est valable pour {= « : 


1 gradn  Qf+ 0: Ë & &) 


n sad 2 +: 2RT? KE) 


En introduisant cette expression dans (77) on 


* trouve 


grade = grad ++ grad®o 


UX 
= ee m(Qr + Q re ©) 


0. re (82) 
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et par conséquent 


grado,, =— (83) 


HU : ; 
3e p(Qi—0@3)grad7. 


L'exemple de la figure r a été calculé à l’aide 
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NOUVEAUX MODES DE FISSION DE L'URANIUM 
TRIPARTITION ET QUADRIPARTITION (suite) 


Par TSIEN SAN-TSIANG, HO ZAH-WEI, R. CHASTEL et L. VIGNERON. 
Laboratoire de Chimie nucléaire du Collège de France. 


V. COMPARAISON DES RÉSULTATS EXPÉRIMENTAUX 
DES DIFFÉRENTS AUTEURS. 


En même temps que les premières publications 
faites par nous (T. H. C. V.) en 1946 sur les preuves 
expérimentales de la tripartition [15] et de la quadri- 
partition [16[, Demers (D) a publié l'observation 
de trajectoires de fission de l’uranium dans l’émulsion 
photographique [17]. 

Il à utilisé une méthode très ingénieuse en dispo- 
Sant une couche mince d’uranate d’ammonium entre 
deux couches d’émulsion, le tout formant une sorte 
de « sandwich » dont les trois couches étaient adhé- 
rentes entre elles. Le point de fission se trouve dans 
la couché d’uranate d’ammonium et peut être 
repéré à 0,7 près (ou o,1 cm d’oir équivalent). 

Il a ne de nombreuses bipartitions et aussi 
six cas dans lesquels une particule légère, qu'il 
suppose être une particule «, Pare d’un point confondu 
avec le point de fission. 

En mars 1947, Green et Livesey (G. L.) utilisant 
la même technique que nous ont publié leurs résul- 


tats sur l’émission d’une particule ISETE lors de 
la fission de l’urarium [18]. Ç 

Dans le même mois, Farwell, Segré et Wiegand 
(E.S. W.) ont conclu par la méthode des coïnci- 
dences, que la troisième particule légère accom- 
pagnant la fission est une particule & [ro]. 


a. Fréquence d'émission du troisième frag- … 
É nombre de tripartitions 
ment léger. — Pour le fHPpOTt 4 
nombre de bipartitions 
les résultats sont les suivants : Es 2 


G., F,S. We. 
0,0029 -E0,0003 0,004 


THE D. 


0,003 E 0,001 0,00/ E 0,002 
; ; ) | ; 


Dans la limite des erreurs expérimentales, notre 
valeur s'accorde avec celles des différents auteurs. 


b. Parcours du troisième fragment léger. — 
Green et Livesey donnent un parcours maximum 


vs mx rte 


PET 
EE le 


AR 


es : 
+ 


tributions statistiques (fig. 15), on voit que 
les courbes ont la même allure, le parcours le plus 
probable indiqué par eux est 27 em alors que nous 
mêmes (T. H. C. V.) donnons 29 em (écart inférieur 
aux erreurs expérimentales). Cependant nous avons 
observé pour les parcours longs (R > 30 cm) une 
fréquence plus grande que G.L. tandis que pour 
les parcours faibles (R < 20 cm) nous avons observé 
* une fréquence plus faible. 


W {en unité arbitraire) 


Fig. 15. — Comparaison de la distribution 
du. parcours de M, obtenue par différents auteurs. 


Demers a indiqué que parmi les six trajectoires 
-qu'il a observées, une trajectoire a un parcours 
Supérieur à 28cm, ce qui est en accord avec nos 
résultats. 

Par”contre, Farwell, Segrè et Wiegand ont donné 
_23 cm comme valeur du parcours maximum, cette 
valeur semble être sous-estimée. 


“ 


c. Angle d'émission du troisième fragment 
léger. —— Le troisième fragment léger est émis 
de préférence dans une direction perpendiculaire 
à celle des deux fragments lourds. Nous avons 
donné (fig. 7b) la distribution statistique des 
angles «, et «, que M, fait avec M, et M, et constaté 
que la courbe de distribution statistique de 
est plus aigüe que celle de l’angle observé avec les 
trajectoires très courtes liées au centre de la trace 
de ‘fissions (phénomènes non calculables qui pour- 
raient être des tripartitions) décrites dans la sec- 
tion (III-B.-h.) (voir la figure 11 b). Contrairement 
à nous, G. L. n’ont pas trouvé de différence notable 
entre ces deux cas, D. et F.S. W. n’ont pas étudié 
éerproblèmez -- 


._ d. Temps d'émission du troisième fragment. 
— La question qui se pose est la suivante : y a-t-il 
un intervalle de temps appréciable entre les deux 
phénomènes suivants 1° fission du noyau d'uranium 
-29 émission du troisième fragment ? : 

F.$S. W. utilisant la méthode des coïncidences, 
ont conclu que cet intervalle de temps est inférieur 
MD. 100$. 
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d avec le nôtre (45 cm). Si l’on compare - Demers a pu indiquer une limite plus serrée 


au plus 2.10 148, 

L'un de nous [29], considérant l'existence d’une 
direction préférenlielle dans l'émission du fragment 
léger, en a conclu que lintervalle de temps, s'il 
existe, est de l’ordre de la vie moyenne du noyau 
composé”; ,;U, Dans le langage courant, on peut donc 
parler de simullanéité de l'émission des trois fragments. 

Feather [30], par une discussion théorique basée 
sur des considérations analogues, conclut que 
l'intervalle de temps est très inférieur à 10745, 

Il semble qu’on soit maintenant généralement 
d'accord avec nos conclusions suivant lesquelles 
les trois fragments sont émis dans un intervalle 
de temps comparable avec la vie moyenne du noyau 
composé. Selon la préférence des auteurs il reste 


_encore une nuance entre l’hypothèse de la tripartilion 


et celle d’une émission secondaire; nous reviendrons 
sur ce point plus loin. 


e. Masse du troisième fragment. — Farwell, 
Segrè et Wiegand ont, d’après des considérations 
d’ionisation, conclu que le troisième fragment est 
une particule «. Leur calcul est basé sur un parcours 
maximum de 23 em d’air pour le troisième fragment; 
étant donné qu’on sait maintenant que ce parcours 
peut atteindre 45 cm (soit à peu près le double) 
et que sa valeur la plus probable est environ 28 cm, 
il semble que leur calcul ait lieu d’être modifié. 

Green et Livesey n’ont pas déterminé les masses 
dans leur publication. D’après une communication 
privée, ils ont calculé les masses tout récemment 
et obtenu pour M, une valeur moyenne d’environ 3,6. 

D'après nos premiers résultats publiés, ainsi que 
par les cas que nous avons vu depuis, il semble que M, 
varie autour de la valeur 5, toutes les valeurs calcu- 
lées sauf une étant compatibles avec une valeur 
unique de M, — 4, 5, 6 ou 7 (voir la note p. 211). Un 
cas spécial où M, — 32 montre que des valeurs très 
éloignées de { sont possibles pour W.. 

En conclusion, t! semble que le troisième fragment 
ait généralement des masses très légères (M; < 10); 
très probablement une partie des cas correspond à 
des particules «. Les résultats expérimentaux ne 
permettent pas d’exclure la possibilité d’avoir 
d’autres noyaux légersque; He. Une très pelile propor- 
lion des cas correspond à des masses M, plus lourdes, 
complètement différentes de celle de la particule «. 
En conséquence, il n’est pas possible de conclure 
si M; a une masse unique, inférieure à 10 dans 
presque tous les cas, constituant ainsi une même 
classe de phénomènes (tandis que les cas, bien moins 
nombreux, où M, est plus lourd, constitueraient 
une classe de phénomènes distincte de la précédente) 
ou bien s’il y a un spectre de masses. 

Il est donc souhaïtable qu'une spectrographie 
de masse des produits de fission soit faite dans la 
région des poids atomiques légers (M < 70) en 
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employant une source intense de neutrons, par 
exemple une pile à uranium: ainsi la question 
pourrait être tranchée. 


Î. Isotope responsable de la tripartition avec 
troisième fragment de long parcours. —— Nous 
avons conclu que *;;U est responsable de l'émission 
du troisième fragment de long parcours. Cette 
conclusion est appuyée par le fait que F.$. W. qui 
utilisaient un échantillon d'uranium enrichi en 
isotope “;,U ont observé la même fréquence du 
phénomène ‘que nous. Ils ont étudié aussi la fission 
de . **’Pu. par des neutrons lents ét observé 


VI. CONSIDÉRATIONS THÉORIQUES. 


À, La HAROrNene 

a. Le modèle de la goutte liquide et la radar — 
D'après le modèle du noyau imaginé par Bohr, com- 
parant le noyau à une goutte liquide, Present [14] a 
signalé le premier qu’une fission en trois fragments 
lourds chargés est dynamiquement possible par la 
capture d’un neutron. L'énergie d’activation serait 
la même que pour la bipartition et l’énergie libérée 
serait d'environ 10 à 20 MeV supérieure à celle de la 
bipartition. 

Considérons une petite déformetion arbitraire de 
la goutte liquide représentant le noyau composé “a Al 
telle que la distance du centre de la goutte à un 
point quelconque de colatitude 0 sur la surface 
soit changée de sa valeur R, à la valeur 


r 


R(B)= Rol1+ao+ a P;(cosû) 
Er + az P;(cos 0) + x; P;(cos0) +...] 


où les «, sont de petites quantités. 

Bohr et Wheeler [13] ont montré que quand 
l'amplitude « du deuxième harmonique de défor- 
mation est importante, il y a possibilité de fission 
en deux fragments à peu près égaux. 

Suivant les mêmes considérations, Present à 
indiqué qu’à une grande amplitude &, du quatrième 
harmonique de déformation correspond un double 
étranglement le long de l’axe de révolution de la 
goutte allongée. Ces deux étranglements vont 
s’accentuer quand la goutte s’allonge et finalement 
il se forme trois gouttelettes alignées. Les grandeurs 
relatives des trois fragments dépendant des ampli- 
tudes et des phases des deuxième et quatrième 
harmoniques de déformation. D’après lui, parfois 
la dissymétrie probable pourrait provoquer la 
séparation d’un des fragments extérieurs en laissant 
les deux autres fragments dans un état excité, 
puis ceux-ci se sépareraient ensuite spontanément. 


b. Interprétation de la légèreté de la masse et de 
la direction d'émission du troisième fragment de long 


une fréquence moitié. onde en phénomène. | 

Dans un travail fait en collaboration avec. 
Mme Faraggi, l’un de nous [27, 31] a montré que la 
fission de *??Th par la capture d’un neutron rapide 
ne donne pas ce genre de phénomènes, résultat 
identique à celui de * 5 SU soumis à l’action des neutrons 
rapides. = 

Tous ces lails nous conduisent à conclure que - 
l'émission d'un troisième fragment de long parcours 
est un phénomène qui ne se produit pas avec tous 
les noyaux susceptibles de subir la fission et dont 
la proportion par rapport à la bipartition n'est pas la. 
même pour lous les éléments qui y donnent lieu. 


j 
| 


parcours. — En utilisant l’image de la goutte: 
liquide et l'hypothèse de la tripartition, l’un de nous 
a pu expliquer la direction de l'émission du troisième 
fragment léger [9]. Nous allons reprendre cette 


Fig, 16. — Schémas simplifiés de la tripartition. Les rayons Î 
des novaux M, et 4: sont calculés d’après la relation. 1 


approchée R =r, X Mi Où ra = 1,5X10-em. Le rayons 
de M, est supposé celui d’une particule œ (R, 9,2 XTOS A cm). 


€, trois noyaux en contac!; 4, M, est en contact avec M 
et M; mais ces derniers sont séparés de 3X r0-1+Cm;. 
e, M, s’écarte très peu de l’axe de révolution du noyau | 
alongé, 


; É 

xpérimentalement, la fréquence de la bipartition * 
est très grande devant celle de la tripartition et 
celle de la quadripartition, il en résulte que l'ampli-# 
tude du deuxième harmonique est en général. très # 


grande par A DADE à carie des monqtes d'ordre 
supérieur. Il semble naturel d'admettre que quand 
l'amplitude du quatrième harmonique atteint une 
intensité suflisante pour provoquer la formation 


d'une gouttelette centrale, cette amplitude soit 
petite et par suite que la gouttelette centrale soit 
de volume petit devant ceux des deux autres gout- 
telettes. 

Supposons que par un effet de vibration violent, 
ou un effet de dissymétrie, M, se déplace légè- 
rement hors de l’axe (fig. 16 b) tandis que les étran- 
glements s’accentuent, on arrive au stade où les 
trois noyaux sont en contact, formant des noyaux 

_ indépendants (fig. 16 c, d). À partir de ce moment, 
les trois noyaux vont être projetés violemment ; 
un des deux noyaux M,, M, (M, pour fixer les 
_idées), pourrait se détacher des deux autres suivant 
le schéma de Present tandis que M, et M, resteraient 
adhérents entre eux et ne se sépareraient qu’un 
peu plus tard. Si l'intervalle de temps séparant 
ces deux phases du phénomène était grand devant 
la vie moyenne, on ne s’expliquerait pas pourquoi 
M, est émis dans une direction préférentielle perpen- 
diculaire à la direction de M.. 

Si au contraire, le temps qui s'écoule est tellement 
court qu'on peut considérer que les trois fragments 
sont séparés presque simultanément, M, va être pro- 
jeté suivant la résultante des forces répulsives F,_; 
et F,_, qui sont proportionnelles aux charges Z, et Z 
des fragments M, et M, et inversement proportion- 
nelles aux carrés des rayons R, et R,. En première 
approximation, nous pouvons prendre en moyenne [/|] 


Lr—90 et Po 0 


eE 


Fu — Vo LOS CN Tel Fo 7;0L10 1 CM. 
Nous sommes ainsi conduits à A0 — 50 (fig. 
ét AO — 140 (fig. 


-rimentale de A0 — 


16 C) 
16 d). La valeur moyenne expé- 
= A — 80 (fig. 16 a) (voir 
fig. 7 à, section III, B, c) s’encadre parfaitement 
entre ces deux limites; cette comparaison semble 
appuyer l'hypothèse de la simultanéité de la sépa- 
ration des trois fragments. 


c. Interprétation de la direction d'émission de la 
particule légère de faible énergie. — Parmi les 
cas de fission de ?5U et *SU (en mélange isotopique 
naturel) par la capture d’un neutron, nous. avons 
observé, ainsi que Green et Livesey, qu'environ 
x pour 100 des traces de fission sont, dans leur partie 
centrale rectiligne, liées avec une trajectoire courte 
de particule légère. La technique photographique 
ne permet pas de déterminer la masse de ces parti- 
cules. Probablement une partie de ces phénomènes 
correspond à la projection d’un noyau de la matière 
de l’émulsion par un fragment de fission juste au 
début de son parcours, phénomène déjà observé 


lette centrale, 
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par F. Toliot” au cours de ses travaux sur Fes rayons 


de recul & [32]. Mais les travaux non publiés de 
Feather et de ses collaborateurs [18] effectués par 
la méthode de coïncidence ainsi que l'existence d’une 
direction préférentielle d'émission (section ITF, B, À) 
semblent indiquer qu'au moins certaines de ces 
trajectoires, surtout celles qui sont émises dans la 
direction perpendiculaire, sont liées au phénomène 
de fission. Le même genre de phénomènes a été 
observé également par l’un de nous en collaboration 
avec Mme Faragoïi dans le cas de la fission du thorium 
provoquée par des neutrons rapides [31]. 

Si nous admettons que ces particules sont émises 
en même temps que deux fragments lourds, l’angle 
moyen observé A+ 260 (fig. 11 b et 14 b) corres- 
pondrait au cas où, au moment de la séparation, 
la gouttelette centrale est très peu écartée de l’axe 
de révolution de la goutte allongée (fig. 16e). 
Les énergies observées de ces particules (1 ou 2 MeV) 
sont très inférieures à celle du troisième fragment 
de long parcours, elles s’accordent qualitativement 
avec le fait que la résultante des forces prévue par 
le schéma simplifié de la tripartition est plus faible 
dans le cas de la figure 16 e que dans le cas de la 
figure 16 €. 


d. Diseussions. — Dans l'interprétation de la 
tripartition par le modèle de la goutte liquide, 
l’amplitude du quatrième harmonique de défor- 
mation peut varier d’un cas à l’autre : il y corres- 
pondrait une variation de grandeur de la goutte- 
c’est-à-dire qu’il semble naturel 
de prévoir un spectre de masses pour M, les petites 
masses étant plus avantagées que les grandes masses. 
Nos observations concordent avec ces prévisions. 

Si nous regardons l’isotope responsable de la 
tripartition, l’accord ne paraît pas aussi satisfaisant : 
a priort, il ne devrait pas y avoir de différence entre 
Nombre de tripartitions 
Nombre de bipartitions 
noyaux étudiés. Les observations expérimentales 
ont montré que l'émission de la particule légère 
de faible énergie est un phénomène commun à 
toutes les fissions étudiées, avec la même fréquence 
(x pour 100 du nombre de bipartitions pour ?#U, 
238U et %2Th); par contre la tripartition avec émission 
d’un troisième fragment de long parcours a seulement 
lieu avec 5Ù et “29pu avec “des fréquences diffé- 
rentes dans ces deux cas. 

D'après le mécanisme de la tripartition déduit 
du modèle de la goutte liquide, nous ne voyons pas 


le rapport » pour les différents 


la raison pour laquelle l’émission d’une particule 


de long parcours n’a lieu qu'avec les noyaux faci- 
lement fissionables 5U et 23?Pu. 

Peut-être celà est-il lié à l’excitation des niveaux 
du noyau composé au-dessus du seuil de fission. 
Il semble que les niveaux moins élevés favoriseraient 
la tripartition avec émission d'une particule de long 
parcours. Cette idée est appuyée sur les faits suivants : 


sont à envisager : 
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19 9Pu est plus facilement fssionable que 285Ù, 
après capture d’un neutron lent (d'énergie ciné- 
tique presque nulle) l'énergie d’excitation comptée 
au-dessus du seuil de la fission sera plus grande 
pour le noyau composé de plutonium que pour le 
noyau composé d'uranium; on constate que la 
fréquence du phénomène est plus grande avec 
l'uranium Œu : 0,002; U : 0,004; nombres obtenus 
avec le même dispositif expérimental par Farwell, 
Segrè et Wiegand). 


20 28U et Th soumis par nous à l’action des 
neutrons rapides produits par les réactions Be-D 
(En = 6,7 MeV) n’ont pas donné de fréquences obser- 
vables du phénomène; l’état d’excitation atteint 
pour le noyau composé, compté au-dessus du seuil 
de fission, était généralement plus élevé ue dans 
les cas cités à l’alinéa précédent. 


Il serait souhaitable, pour vérifier notre suppo- 
sition, que des expériences soient entreprises avec 
des isotopes séparés. #SU et PU. dinsi qu'avec #°Pu 
soumis à des neutrons monocinéliques de différentes 
énergies. 


B. L'émission d’une particule «. 


a. Réaction (n, «), tripartition, ou radioactivité x ? — 
L'aspect de l’ionisation du troisième fragment 
léger dé long parcours semble indiquer qu’il s’agit 
d’une particule «, bien connue dans les réactions 
nucléaires. Les expériences de Farwell, Segrè et 
Wiegand appuient cette hypothèse et nos obser- 
vations, dans la majorité des cas (M, < 10), sont 
aussi compatibles avec une masse unique égale à / 
(section III, B, d). 

L’étude de la direction d'émission de M, permet 
de donner des renseignements sur la nature intime 
du phénomène. 

S'il s’agit d’une particule «, trois possibilités 
10 réaction (n, «) [33]; 2° tripar- 
tition [29]; 3° radioactivité x [30, 33], suivant que 
la particule « est émise, avant, pendant ou après 


: Ja fission. 


Si elle était émise avant la division du noyau (1°), 
sa direction serait de préférence suivant l’axe de 
la goutte déformée, car la barrière de potentiel est 
moins élevée dans cette direction [29]. 

Si elle était émise, après la fission (3°), par un des 
fragments lourds, sa direction aurait une distribution 
au hasard par rapport à ce fragment lourd [20]. 


VII. CONCLUSIONS. 


Utilisant la technique de la plaque photographique, 
nous avons mis en évidence la tripartition et la 
quadripartition du noyau composé ?U,. 

Dans le cas. de la tripartition, les deux fragments 
lourds ont des parcours et des masses semblables 
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Ja fission, tandis que si on parle de radioactivité « 


La direction précrenttelles “observée pour M est; 


perpendiculaire aux directions de M, et de M, celà K 
indique que si la particule & part avant ou après : 
la fission des deux fragments lourds, l'intervalle: 


de temps séparant les deux phénomènes est vrai- 


semblablement très court[29l; il n'y a par conséquent 


pas possibilité de distinguer entre une réaction 


du type (n, à) et une radioactivité « d’un des frag- 
ments de fission. ; 
Tout ce que nous pouvons déduire de l’expérience, 
c'est que la particule à est certainement émise sous 
l'influence presque égale et simultanée des champs 
coulombiens des deux noyaux lourds M, cet M, 


l'émission d’une particule « ne constituerait qu'un cas 


particulier de tripartition. 
La seule nuance entre l'hypothèse de la tripar- 
tition et celle d’une radioactivité « serait que si 


l’on parle de tripartition, on considère la goutte- 


lette centrale comme préformée au moment de 


(ou émission secondaire) on considère que la parti- 
cule « se forme aux dépens d’un morceau de fission 
immédiatement après la séparation des deux mor- 
ceaux lourds, cette dernière façon de voir s'accorde 
avec celle de Present dans son explication de la 
tripartition. Nous pensons qu ‘entre les deux façons 
de voir il n'y a qu'une différence de phraséologie 
puisqu'on ne peut pas préciser le temps très court 
qui séparerait les deux phénomènes successifs. 


b. Discussion sur une suggestion de Feather. — 
En considérant la stabilité des noyaux excités 
vis-à-vis de la radioactivité «, Feather [30] a suggéré 
que l’émission des trajectoires courtes (de faible 


énergie) serait due aux produits de fission de Z æ 40 


tandis que l'émission des trajectoires longues (de 
grande énergie) serait due aux produits de fission 
de Z = 60. Cette su ggestion paraît très séduisante 


par sa simplicité, mais les faits suivants paraissent E 


la contredire : È 


10 La masse moyenne du fragment M, de tripar- 
tition est 130, d'après l’étude des produits de. 
fission de ##U par les neutrons lents [4], ces masses 
correspondent plutôt à Z &w 5o qu'à Z = 60. Les 
masses correspondantes de Z  6o étant environ 150. 

20 La masse du cas particulier M, — 
tout à fait hors de la zone M, == 4. - 


Il semble qu’au moins une partie des faits expé- 


rimentaux ne-peut pas être expliquée pu la sugges- 


tion de Keather. 


à ceux qu'on observe avec la bipartition. Le troi- 


sième fragment a le plus souvent une masse infé- 


rieure à 10, mais il est arrivé que cette masse atteigne 
la valeur 3°. Le parcours du troisième fragment 
cm à 45 cm d'air avec la valeur la plus. 


varie de 2 


Ed 


32 este 


; F à CE { 
ET PET EP PAS 


ù ne 


© probable 28 es cm; sa on est de cs 


_ perpendiculaire à celles des deux fragments lourds. 
Les masses M, semblent former un spectre avec - 


un maximum de probabilité vers la masse très 
légère M, = 5, cependant, en tenant compte des 
erreurs expérimentales, tous les cas sauf un sont 


. compatibles avec une valeur unique de la masse 


comprise entre 4 et 7. L'énergie cinétique de tripar- 
tition 154 + 10 MeV, légèrement supérieure à celle 


de la bipartition, est qualitativement en accord 


avec la prévision théorique. La fréquence du phéno- 


mène par rapport à celle de la bipartition est 


0,003 = 0,001. 

Un groupe de particules de faible énergie (énergie 
probable : 1 à 2 MeV) semble être émis au moment 
"de la fission. S'il entre dans la catégorie des phéno- 
mènes de tripartition, sa fréquence serait o,o1 de 
celle de la bipartition. 

Deux cas de quadripartition ont été enregistrés; 
l’un correspond à quatre fragments lourds (M > 20), 


l’autre correspond à trois fragments lourds et un 


léger. 


L'énergie cinétique moyenne est environ 
100 MeV, en bon accord avec la prévision théorique. 


_ La fréquence de la quadripartition par rapport à 
_ celle de la bipartition est 0,0002 + 0,000 15. 


Dans le cas de la fission du noyau composé 2#U, 
ni tripartition avec troisième fragment de long 


Parcours, ni quadripartition n’ont été observées. 


Seule l'émission de la particule légère de parcours 
court a été enregistrée. Son énergie, sa fréquence 
et l'aspect du phénomène sont tout à fait semblables 
à ce qui a été observé avec le noyau composé *#U. 
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F Une partie de ces résultats est confirmée par 
d’autres auteurs. 

La discussion de la plupart des résultats concorde 
d’une façon satisfaisante avec le mécanisme de 
la tripartition déduit du modèle de la goutte liquide 
du noyau (les trois fragments se sont séparés presque 
simultanément). Le phénomène de tripartition avec 


l'émission du troisième fragment de long parcours 


pourrait être lié avec les niveaux moins élevés 
du noyau composé. 
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APPENDICE N° 1. AFFAISSEMENT DE LA GÉLATINE APRÈS LE FIXAGE. 


Les émulsions que nous avons employées ont 
sensiblement la composition suivante (d’après les 
indications du fabriquant Ilford Ltd) : 


S2 pour 100 
15 » 


Après développement et fixage, dans les conditions 
où nous avons opéré, la plaque reste pratiquement 
transparente; une partie infime des grains de bro- 
mure d'argent a été transformée en grains d’argent 
réduit. 

Le départ du bromure d'argent a pour conséquence 
une diminution du volume de la couche de gélatine. 
Cette couche étant adhérente au verre de la plaque, 
la diminution de volume ne peut se faire que dans 
le sens vertical : la gélatine subira donc une sorte 
d’ « affaissement ». 

Il est facile de constater que cet affaissement se 
fait de façon régulière et n’est pas accompagné de 
déformations locales dans le plan horizontal. En 
effet, la trace laissée dans l’émulsion par un rayon « 
du polonium par exemple, reste rectiligne au cours 
des’ Dee de fixage et de développement, et ceci 


qu'elle soit peu ou très inclinée sur l’horizontale : 
on le constate immédiatement en observant une 
plaque au microscope. On constate aussi au micro- 
scope que la surface de la gélatine après fixage et 
développement est restée horizontale. 


a. Reconstitution des figures dans l’espace. —— 
Soit un phénomène nucléaire (choc d’une particule 
en projetant une autre, tripartition, quadriparti- 
tion, etc.). Les longueurs qui ne sont pas horizon- 
tales, les angles qui n’ont pas leurs deux côtés 
horizontaux, subissent une déformation due à l’affais- 
sement de la gélatine, il y aura lieu de reconstituer la 
figure dans l’espace telle qu’elle était dans la gélatine 


avant le développement, on opèrera en s'inspirant 


des règles suivantes : 

x. les projections horizontales des longueurs se 
conservent; 

GB. les projections horizontales des angles se 
conservent; 

y. pour reconstituer la figure, il faut multiplier les 
cotes verticales mesurées après fixage par un fac- 
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A 


teur k supérieur à l’unité et qui nè dépend que de 
l’émulsion utilisée. 


b. Calcul approché du facteur k. — Connais- 
sant la densité du bromure d’argent cristallisé (6,46) 
et celle de la gélatine sèche (1,40 environ), nous 
avons, pour 1 g d’émulsion, : 

, 82 

Volume du bromure d'argent... AE = 0, 125710m? 
», 40 
0,18 


» delarvélatine tee = 20) 120 CE 
1,4 


» der g d’émulsion...... 0,256 cms 


sachant qu’un volume de 0,129 après fixage corres- 
pond à un volume de 0,256 d’émulsion, le facteur k 
par lequel il faut multiplier les cotes verticales 
obtenues est égal à 

0,256 


Ù — = 1,97 
0,129 Ru 


(à 10 pour 100 près). 


Les causes d’erreur sont l'incertitude sur les pro- 
portions exactes de Br Ag et de gélatine et l’incer- 
titude sur la densité de la gélatine sèche suivant sa 
provenance. 

Une couche d’ ation d’une épaisseur de 4op 
avant fixage, verra celui-ci réduire son épaisseur 


UT . . 
à ; X 4op, soit 20p. environ. 


c. Détermination précise des facteurs k' et k 
de multiplication des cotes verticales. — Le 
calcul précédent présente, comme nous l’avons vu, 
une certaine incertitude. Nous avons donc irradié 
dans le vide deux plaques « nuclear research C, » 
avec des rayons « provenant d’une source de polo- 
nium. 

Dans une de ces plaques, les rayons & arrivaient 
tangentiellement à la surface de l’émulsion, une 
statistique des longueurs a donné, pour leur parcours, 
la valeur moyenne suivante : ! — 20,3 divisions du 
micromètre oculaire (à 0,5 pour 100 près). 

Dans la deuxième plaque, les rayons « afrivaient 
perpendiculairement à la surface de l’émulsion. 
Nous avons constaté qu'il était possible de suivre la 
trajectoire verticale au microscope. Leur longueur 
était mesurée au moyen d’un tambour gradué 


APPENDICE N° 2. — 


DISTINCTION D'AVEC LES TRIPARTITIONS. 


Lorsqu'on se trouve en présence de trois traces 


issues d’un même point de la gélatine (fig. 3), 


on peut être en présence d’une tripartition; en ee cas, 
les trois branches jouent des rôles symétriques. 

I ne faut pas se hâter de conclure, car on peut être 
en présence d’un simple choc nucléaire (ici les trois 
branches ne jouent plus des rôles semblables). 
La figure 17 représente un tel cas. 


actionnant le mouvement Vertical lent du microscope. 
Une statistique a donné, pour la longueur des 
traces, la valeur suivante : l' — 4,88 division de la 
roue (à 2 pour 100 près). 

Il en résulte que si une cote verticale après fivage 
a une mesure de x divisions de la roue du mouvement 
lent, cette cote avant le développement avait la valeur 
suivante : k! x divisions du micromètre oculaire, avec. 


. 


20; 


k = 
(è Ts 


A émis SRE Amir rt 
acte Le-tinot h m ni ue LAS 


= 4,16 (à 3 pour 100  prèsi. 


œ 


2 


C’est sous cette forme que le résultat est le plus 
utilisable. On peut se proposer de calculer le facteur k 
défini dans le paragraphe précédent. Une division 
du micromètre oculaire représente 1,05u. La vis 
micrométrique avait été soigneusement étalonnée et 
au point où nous l’utilisions : une division de la 
roue vaut 1,87. Nous avions donc 


CPR 7 


CESR PE DES PESNP DUT 


} 
20,3 X 1,09 L—KXx<4,88 X1,87-U, 
d’où 

k—2,33 à 10 pour 100 près. 


L'erreur relative sur k est supérieure à l’erreur sur k' 
car l’on ne connaît pas avec une grande précision 
la valeur en microns des divisions du micromètre 4 
oculaire ou de la roue verticale du mouvement lent. 

On constate que les valeurs obtenues pour le à 
facteur k par calcul approché et par mesure directe 
de longueurs de traces sont compatibles. 

Nos mesures et nos statistiques dans le plan 
horizontal sont effectuées avec l'unité suivante : 
division du micromètre oculaire. Ce n'est qu’ensuite 
que nous effectuons une transformation en microns. 
De même, l'unité de déplacement vertical est la 
division de la roue du mouvement lent, le micron, 
là aussi, n’est qu’une unité secondaire, il était done 
plus commode pour nous d'utiliser le facteur k' 
mesuré avec une meilleure précision. Mais ce facteur 
n'est valable que pour le microscope que nous avons 
utilisé. Au contraire, le facteur k est valable d’une 
façon générale quel que soit le microscope employé. 

Des mesures statistiques effectuées ‘avec des 
rayons « dont la pente par rapport à l'horizontale 
était 45 ou 60° ont donné les mêmes valeurs pour k! 
que les mesures avec une incidence verticale. - 


CHOCS NUCLÉAIRES. 


dés haréarues dsl hé rritiare apte at teinté 


irteeteiimersehnhste 


L 


Un atome d'uranium s’est seindé au point (0) 
en deux fragments de masses M et N. Le fragment N 
s'est arrêté en T. Le fragment M, lancé en sens 
inverse en direction OR, a rencontré un noyau 
cn P; dévié de son chemin, il s’est arrêté en S. 
Un noyau de masse MX primitivement au repos 
en P a été projeté par le choc, a laissé la trace PO 
et s’est arrêté en Q. 


eo 


adm p X dess 
Pare rad tons 


© NOUVEAUX MODES DE FISS SION DE L'URANIUM 


Il est. den que rien ne signale le point O où 
s’est produit la fission, il est « quelque part entre T 
Let Pr 


. Fig. 17. — En O : fssion d’un noyau U en deux masses M 
et N. En P : choc de M contre un noyau de masse M* 
précédemment immobile, 


Appelons : 


_ V , la vitesse de M au point P, juste avant le choc; 
E, l'énergie cinétique correspondante ; 
V', la vitesse de M au point P, juste après le choc; 
E", l'énergie cinétique correspondante; 
V*, la vitesse de M* au point P, juste après le choc; 


3 EY* l'énergie cinétique correspondante; 
: M 
P; le rapport des masses p — DLÉ 


g, le rapport des vitesses des deux particules 


TL. 


juste après le choc g — 2 

_«, B, y, les angles (voir fig. 17). 
a. Application des principes de la méca- 
nique. — L'application au point P des principes 


de conservation de l'énergie cinétique et de l'impulsion 
conduit facilement aux résultats suivants : 


1° Relation entre les masses 
M sin? y — sin? f 
sin? & 


29 Relations entre les vitesses 


[4 sin y Sinœ 
ps; PERRE sin? 6? 
v” . sinÿ sinc Herr sin 
Pas EE ; rt 
ARE 2: sin y sin2f? Prin e 


3° Relations entre les énergies cinétiques 


E sin? y 
SE sin? y — sin? 6” 
1e sin? 5 pe sinae 
Er © sin? y — sin G° ET sin?f 


Il est très remarquable que fous ces rapports 
puissent se déduire uniquement de la connaissance 
des qe de la figure (sous réserve qu'on sache quelle 
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branche TOP correspond à la particule incidente 
avant choc, quelle branche PS correspond à la 
particule déviée et quelle branche PQ correspond 
au noyau M* projeté). 


4° L’angle B est toujours aigu; 
59 L’angle + est toujours supérieur à l'angle B;, 
60 Si M* < M, l'angle à« + B des deux particules 


après le choc est toujours aigu. Si M* — M, 
l’angle 4 + 6 est toujours égal à . St M*>.M, 


l'angle & +8 est toujours obtus (£ <a +p< r). 


b. Étude des différents chocs possibles. 
l’on se trouve en présence d’une étoile à trois 
branches, dans le cas général on ne sait pas quel rôle 
attribuer à chacune (il y a trois rôles possibles pour 
une branche : être le segment TP, le segment PS 
ou le segment PQ). Finalement, six figures de choc 
différentes sont possibles. 


Fig, 18. — La particule projetante ne peut se déplacer 
dans le sens indiqué par les flèches de la figure 18 a car la 
condition y > $ n’est pas vérifiée. Le sens de parcours 
indiqué par les flèches de la figure 18 b vérifie y > £. 


On démontre facilement que la condition angu- 
laire 7 <8B élimine dans tous les cas trois de ces 
figures : ilreste donc trois figures de choc possibles. 


Simplifications : Parfois on a des raisons sûres 
de connaître le rôle d’un segment; par exemple, 
si les particules M proviennent toutes d’une même 
source, on reconnaîtra la particule ‘incidente à sa 
direction. Dans l'étude de nos plaques, souvent 
deux branches de l'étoile étaient épaisses et la 
troisième avait une ionisation bien plus faible : elle 
ne pouvait pas être la trace d’un fragment de bipar- 
tition, c'était donc à coup sûr la particule projetée 
(si l'apparence était due à un choc); on se trouvait 
alors en présence de deux figures de choc possible, 
mais on démontre facilement que la condition angu- 
laire + = 8 en élimine alors une dans tous les cas 
et qu'il ne reste qu'une figure de choc à étudier 


(fig. 18). 


c. Élimination du cas du choc. - Lorsqu'on 
sait quelle est la particule incidente, quelle est la 
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particule émergente et quel est le noyat Hot, 


les angles caractéristiques x, B, y peuvent être lus 
sur la figure (fig. 18 b). 


10 Les rapports pet q peuvent être déterminés immé- 
diatement d’après les mesures de &, 5, y par les 
formules du paragraphe a. 


20 Choix provisoire d’une masse M*. — M* ne 
peut avoir qu’une des valeurs suivantes : 1, 12, 14, 
16, 80, 108, 235-238 (masses des noyaux H, C, N, 
O, Br, Ag, U). Dans le cas (fréquent dans notre 
travail), où l’ionisation due au passage de M* est 
nettement plus faible que celle d’un fragment de 
fission et nettement plus forte que celle d’un proton, 
il ne reste que les valeurs 12, 14, 16. Si, au contraire, 
M* avait produit une trace épaisse, dans le genre 
de celles des fragments de fission, on ne pourrait 


prendre pour valeurs possibles de M* que 80, 108, 238. 


On va alors choisir provisoirement une masse M*, 
parmi celles qui sont possibles, par exemple M* =:59 
et l’on va rechercher si la figure observée peut corres- 
pondre à La projection d’un noyau de carbone par un 
fragment de fission. 


59 Caleul de M par la formule M—pM*. — 
Ce calcul est immédiat, puisqu'on a déterminé p 
au paragraphe 1°. Supposons que p = 20, M* — 1», 
on aurait M — 20 X 12 — 240, masse absolument 
excessive pour être celle d’un fragment de bipar- 
tition. Si.p —2, on trouverait M —2 Xu2= 514, 
valeur absolument insuffisante pour être la masse 
d’un fragment de bipartition. Dans ces deux cas, 


on pourrait, sans aller plus loin, éliminer le cas d'un 


choc contre un noyau de carbone. 

Il arrive qu’on obtienne pour M des valeurs 
tout à fait vraisemblables; ainsi, si p calculé au 
paragraphe 1° était égal à 10, la masse M — pM* 
serait 120, valeur admissible. 


4° Confrontation des vitesses V' obtenues de 
deux façons différentes. — Si le calcul du 
paragraphe 30 donne pour M une valeur vraisem- 
blable, il faut pousser plus loin la discussion. Connais- 
sant la valeur provisoirement admise pour M*, 
et connaissant le parcours restant (longueur de la 
trace PR de M*) on déduit immédiatement V*. 
Le rapport q a été calculé au paragraphe 1°, on 
calcule une première valeur de V' par la for- 
mule V'=qV*, 


APPENDICE N° 3. -- CALCUL DES TRIPARTITIONS ET DES QUADRIPARTITIONS 
PAR APPROXIMATIONS SUCCESSIVES. 


Reportons-nous au paragraphe a de la section ITI-B 
et à la figure 3 représentant une tripartition; on s’est 
donné arbitrairement deux axec rectangulaires 0%, Oy 
dans cette figure. 


rence expérimentale, on ne peut alors rien conclure 


“D'un äutre côté, on a dun au paragraphe 30 
la valeur de M, connaissant d’autre part la PE 
de la trace PS de cette particule après le choc, . 
on obtient immédiatement une deuxième valeur de V!. « 

Si les deux valeurs obtenues concordent dans Ja 
limite des erreurs expérimentales possibles, la A 
figure peut être due à un choc. 4 


50 Essais de nouvelles valeurs pour M*. —. 
Il peut arriver qu’en appliquant les opérations des . 
paragraphes 30 et 4° à une masse provisoirement . 
choisie, on arrive dans tous les cas à une absurdité. - 
La valeur choisie pour M*, 12 par exemple, ne peut . 
pas convenir, mais il n’est pas dit qu’une autre « 
valeur, 16 par exemple (noyau d’ oxygène), ne puisse 
pas convenir. On recommencera donc les opérations … 
des paragraphes 3° et 49 avec MX — 16. TE 


60 Il faut aussi tenir compte du fait que les 
rapports p et q ne sont pas connus avec une très. 
grande précision puisqu'ils dépendent des angles lus « 
sur la figure. Si l’angle « est petit, de faibles varia- « 
tions absolues de « causeront d’énormes “ariations 
sur p et q et il faudra en tenir compte en reprenant 
la discussion pour différentes valeurs admissibles « 
pour « (par exemple faire la discussion pour & — 29, 39, M 
4°, 59, 6°, 79 si la valeur mesurée est x — 40,5 + 20 5) 


7° Conclusion de la diseussion. — a. Si ape 
la discussion menée comme il vient d'être dit, le | 
choc n’est possible dans aucun cas, on a donc affaire 
à une tripartition. b. Dans certains cas, d’après la. 
discussion ci-dessus, un choc peut interpréter l’appa- 


cela peut être un choc, cela peut être une tripar-” 
tition. "a 

Il y aurait alors lieu de calculer N — »36— M, 
masse du deuxième fragment de bipartition, et de. 
voir si là branche PT correspondante a un parcours 
trop court (ou bien très excessif); on arriverait. 
encore à éliminer quelques cas mais on conçoit que 
certaines apparences pourraient encore s’inter- « 
préter de deux façons, sans qu’on puisse conclure 
entre le choc et la tripartition. Quand de tel cas se 
sont présentés, nous n’en avons pas tenu compte M 
dans la statistique des tripartitions. “à 

Dans tous les cas où la trace légère avait un 
parcours supérieur à {Cm d’air équivalent, nous 
avons constaté que l'hypothèse du choc conduisait 
à une incompatibilité. Il y avait donc tripartition. 


Pour les trois valeurs R,, R,, R, (longueurs en 
centimètres d'air équivalent des trois traces), nous! 
construisons les trois courbes donnant la vitesse en“ 
fonction de la masse ainsi qu’il a été dit en ILI-B a 


_ successives. 


SNS 


a. * Systèree de trois équations donnant les M. 
— Considérons les conditions de masse, ZM — 236, 


“et d’ impulsion, 2M (ie — 0, définies en III-B a. 
Cette dernière équation est une équation vecto- 

rielle; par projection sur les deux axes Ox, Oy, 

on voit qu'elle est équivalente aux deux premières 


équations du système ci-dessous : 


Mio Vie MB Pr My Vs = Pro, 
M, C4 " + M 6’ la At M Y' Vs P; = 0,, 
M, + Mo + M; 236: 


Système A 


% a; 6, B'; y, y’ sont les cosinus directeurs des direc- 
tions des trois traces. Dans la figure 3, on a les 
valeurs numériques suivantes : 


QE IA NP — 
B=— sind ——0,174, {= 


sin 0; = 0,960, 
B'= cos 0 — 0,985, Y— 00503 —10,25. 

Si l'on considère les V comme connus, on remarque 
que le système À comprend trois équations linéaires 
à trois inconnues M,, M,, M. On peut le résoudre 
par rapport à ces inconnues. 

Cette remarque nous a inspiré 
d’approximations successives suivante : 


la méthode 


b. Masses et vitesses par approximations 
10 Départ. — Commençons par 
prendre trois valeurs arbitraires pour les masses M.,, 
M;, M;. Afin de réduire le nombre des approxima- 
tions, on s’efforcera qu’elles n'aient pas des valeurs 
trop éloignées de la solution et pour cela on se 
laissera guider dans le choix par des considérations 
d’ionisation des trois particules. 


= 


20 Première approximation. — . Vitesses. 
On lit.les vitesses V,, V,, V,; en regard des masses 
sur les courbes V—}M correspondant respecti- 


- vement aux parcours R;, R:, R. 


6. Masses. — Portant ces valeurs approchées des 


vitesses dans le système A, on le résout et l’on obtient 


pour les masses des valeurs mieux approchées que 


celles qu’on avait prises au départ. 


3° Deuxième approximation. — %. Vitesses. 
On lit les vitesses sur les courbes V — M en prenant 


pour valeurs des masses celles qu’on a obtenues à la 


fin de l’approximation précédente. 


6. Masses. — Portant ces vitesses dans le sys- 
tème A, on résout celui-ci et l’on a des valeurs des 


_ masses meilleures que dans l’approximation précé- 


 dente. 


4° Approximations suivantes. Arrêt. — On continue 
les approximations suivant le même principe. On 
constate qu’elles convergent rapidement : par exemple, 


_ les masses de la 42 approximation coïncident, à 


moins d’une unité de masse près avec celles de la 
3e approximation. On peut alors s'arrêter et, en pre- 
nant les entiers les plus voisins des masses obtenues, 
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on a la solution. En effet, le calcul algébrique condui- 
rait, le plus souvent, à des masses non entières. 
On constate que 


me — M, (o4 Vi cu M; B Ve y M; de V; 


qui représente la projection P, de la résultante des 
impulsions sur l’axe Ox est alors inférieure à 
1 pour 100 de l’impulsion M, V, d’une des particules 
lourdes. Il en est de même de la projection de cette 
résultante sur Oy. 

La précision est limitée par le fait que le manque 
de données expérimentales sur les relations V—ÆR 
pour certaines zones de masses entraîne une impré- 
cision sur les courbes V — M que nous avons cons- 
truites. De plus, il y a une certaine fluctuation 
individuelle des parcours pour une même vitesse 
(ou de la vitesse pour un même parcours). 


5° Nous donnons ci-dessous un exemple numérique 
des valeurs obtenues dans le cas du calcul de la 
tripartition de la figure 3 : 


M. M: M, Va. Pa. . 
Valeurs arbi- 
traires prises 
au départ.... 116 116 “ 
1*e approximation : 
Vitésses 74 DE, ONLINE 20 
Masses se 122 104,6 9,4 
22 approximation : 
D'ITESSeS TEE 10,9127% 2950 
Masses:o.: 2: 129,00 P109 104 
3° approæimation : 
Vitesses...:.. TO DUT 002085 
Masses. tent 124101025020 ,3 


4e approæimation 5 


Les valeurs des vitesses sont les mêmrs que dans 
la 3° approximation, il en est par suite de même 
des masses. 


On constate que dans la 17e approximation, la solu- 
tion était atteinte à 2 pour 100 près des valeurs des 
masses. 

Solution admise définitivement : 


Mi = 124; Ms — 703, M3 = 9. 


Valeur de M, V, en unités arbitraires : 1352. 
p« dans les mêmes unités : —9 —0,7 pour 100 


de M, V;; 

.p, dans les mêmes unités : 9 — 0,7 pour 100 

de M, Han 

Module de De — È2M Ÿ dans les mêmes unités : 
MOENOLE 


soit 1 pour 100 de M, V.. 
émis. 


14. 


6° Cas où des neutrons ‘seraient 


Imaginons qu’on connaisse leur nombre ainsi que 
la direction de leur vitesse, le calcul ci-dessus reste- 
rait valable à condition qu’on introduise un deuxième 
membre non nul (projection sur Ox ou Oy de l’impul- 
sion totale des neutrons) dans les deux premières 
équations du système A et qu’on remplace dans 
la troisième équation 236 par 233, 234 ou 9235, 
suivant qu'il y aurait eu 3, 2 ou 1 neutrons émis. 
Il est impossible de connaître ces corrections et l’on 
devra résoudre le système À comme si aucun neutron 
n’était émis. Remarquons qu’en appelant À l’impul- 
sion d’un neutron de 4MeV, l'impulsion d’une 
particule « de 16 MeV sera 4 À, l'impulsion d’un 
fragment de fission ER masse 140 et d’énergie 60 MeV 
sera 46 À. \ 


On en déduit que l'erreur commise en négligeant 
les neutrons qui pourraient être émis introduit, 
dans le calcul, une erreur de 2 pour 100 au plus sur 
la masse des deux fragments lourds de la tripar- 
tition et une erreur qui peut atteindre environ 
une unité de masse sur la masse du fragment léger: 


c. Calcul des quadripartitions par approxi- 
mations successives. — Nous pensons que, dans 
le cas général, il n’y.a aucune raison pour que les 
quatre traces soient coplanaires. On considérera 
trois axes rectangulaires de coordonnées, le sys- 


j 
‘ 


tement équivalente à celle qui a été exposée’ dans 
l'Appendice n° 3; l’emploi de constructions gra- 
phiques évite d’avoir plusieurs fois de suite à 
résoudre algébriquement un système de trois équa- 
tions à trois inconnues. 


On commencera par construire les trois courbes 
donnant l’impulsion en fonction de la masse pour 
les trois parcours R;, R:, R:. 

Les projections p. et p, de l'impulsion totale 


. 
p=>MV. 


sur Ox et Oy sont 


Px=d@ Pit pr+ y ps. 
Py= pat PpatY'pa, 


Pis Pas Ps étant respectivement les impulsions M, V,, 
M: Va M3 Vs à, 6, y, «x, P', y’ étant les cosinus 
directeurs des directions des trois traces par rapport 
aux axes Ox et Oy de la figure 3). 

Si les valeurs des masses M,, M,, M, sont cor- 
rectes, en déterminant les impulsions p;, Pa, Ps au 
moyen des courbes MV —#f(V), on trouvera, 
pour p. et p, des valeurs nulles et, par suite, le 
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.ie nombre des équations du système A resterait 


APPENDICE N° 4. — CALCUL DES TRIPARTITIONS 
PAR UNE MÉTHODE GRAPHIQUE. 


La méthode que nous allons exposer est stric- 


tème À comportera une ten de plus Cprojéction 
de l’impulsion résultante sur Oz — o). 

Il y a une inconnue de plus (M). 

Le système se résoudra en suivant exactement 
une marche calquée sur celle qu’on utilise dans les 
tripartitions. La solution compatible avec l'aspect 
expérimental sera unique. 

Si les quatre particules sont coplanaires, une des 
quatre équations du système A disparaît, il y a donc 
quatre inconnues et trois équations, on pourra se 
donner une valeur arbitraire pour M, et calculer les 
trois autres masses; on obtiendra ainsi une première 
solution. Ensuite, on se donnera d’autres valeurs 
arbitraires pour M,, il y aura une simple infinité de 
solutions compatibles avec l’aspect expérimental. 

En fait, dans la première quadripartition que nous 
avons observée, les quatre traces étaient sensiblement 
coplanaires. Dans la deuxième, les trois traces 
lourdes étaient presque coplanaires, cependant 
l’aspect de l’ionisation (et aussi des considérations 
géométriques) imposait 3 < M, < 6; dans ce cas, une 
solution unique a pu être calculée en posant M, —4, 
l'erreur commise sur M, ne retentissant pas beau- 
coup sur les autres masses. 

S'il y avait plus de quatre particules émises, : 
le nombre des inconnues augmenterait, tandis que 


: 


1 
} 


égal à 4. Il y aurait donc indétermination. 


module de l’impulsion résultante sera nul 


F | = p=W(pz}+ (py} = 0. 


0 75 130 135 140 


45 M; 
Fig. 19, — Y, y’: courbes d’égal module 
de l’impulsion résultante p; 

Z : point où p = o. 


correspondent à la ie P- et p, ne seront pas 
simultanément nuls, p sera différent de zéro. 

Prenons une feuille de papier quadrillé (fig. 19) 
_avec deux axes OM,, OM,. Au point A de la figure 
correspondent M, — 100, M, — 120 ef, par suite, 
M3 = 236 — (M, + M;) — 16 (condition de masse). 

Prenons comme unité des impulsions, celle qu’au- 
rait une particule de masse 1 animée d’une vitesse 
de 10°cm:s. L’impulsion d’un fragment lourd de 
_ bipartition est de l’ordre de 1300 unités, cellé d’une 
particule « de 20 MeV est 122 unités, celle d’un 
neutron de {4 MeV est 27 unités, 

Calculons p pour le point A, soit p — 130 pour 
fixer les idées, inscrivons 130 à côté du point A. 
Opérons de même pour un certain nombre de 
points B, C, D, E, F, G, H, I, etc. du graphique. 
Entre ces points, il est facile de faire des interpo- 
 lations paraboliques graphiques et de déterminer 
des points q, r, s, {, u, v, où p est égal à 5o unités. 
Ces points seront joints par une courbe +. 

__ La courbe y est le lieu des points tels que le 


module de p soit égal à 5o unités. C’est une courbe 
ferrée. La courbe y’, lieu des points où p est égal 
à 150 unités entoure la courbe Y. 

Le point Z correspondant à la solution sera entouré 
par la courbe y, cette dernière ayant une forme régu- 
Hère, on apprécie facilement la position de Z quand 
les courbes y et y’ ont été tracées. 

On calculera, à titre de vérification, la valeur 


at 
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de p pour les masses M,, M,, M, correspondant au 
point Z qu’on vient de déterminer graphiquement, 
on s’estimera satisfait si p < 25, c’est-à-dire si la 
compensation des impulsions a lieu à 2 pour 100 
près de l’impulsion d’un fragment de bipartition. 
Si cela n’a pas lieu, une légère retouche du point Z 
faisant varier les masses de 1 ou 2 unités fournira 
la solution. 


L'avantage de cette méthode graphique est qu’on 
se rend immédiatement compte de l’ordre de gran- 
deur de la précision (quantité dont on peut faire 
varier les masses si l’on admet que p atteigne la 
valeur 25, soit 2 pour 100 de l’impulsion d’un frag- 
ment de bipartition). Sur la figure 19, toute combi- 
naison de masses figurée par un point intérieur à 
la courbe +” (p — 25) peut être considérée comme 
possible. 


Remarque. — Dans une quadripartition, il y a 
une masse M, en plus, cette méthode ne s'applique 
donc pas, car il faudrait faire un graphique OM,M,.M; 
dans l’espace, ce qui ne serait pas commode. 


Note ajoutée à la correction. — Dans le calcul de la masse 
de M;, nous avons utilisé, comme première approximation, 
la relation entre R et V pour la particule +. Étant donnée la 


relation théorique R — BI V) et la charge Z de M, étant 


inconnue, nous avons introduit une erreur incontrôlable sur 
la masse de M,. Pour la valeur calculée de M, 4, l'erreur 
estimée serait environ une unité de masse. 


Manuscrit reçu le ro juillet 1947. 
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DÉTERMINATION DE LA VALEUR LA PLUS PROBABLE DES GRANDEURS STATISTIQUES 
I. GÉNÉRALITÉS 


Par GEorGes ALLARD. 
Sommaire. — Résolution du problème suivant : soit æ une certaine grandeur mesurable pour des 


particules individuelles; la probabilité pour que x soit compris entre x, et x, dépend de [ para- 
., «3; des expériences effectuées sur un grand nombre de particules montrent qu’ilyenam 


mètres œ, @, «. 


pour lesquelles x est compris entre x, et x, n, pour lesquelles x est compris entre x, et æ,, ...; déter- 
miner la valeur la plus probable des { paramètres œ, &, ..., æ et les écarts quadratiques moyens 


correspondant da, ..., da}. 


Dans un mémoire des Proceedings of the Royal 
Society of London (1), Peiïerls a étudié la précision 
des mesures effectuées au moyen de compteurs. 
Il m'a paru intéressant de reprendre ce travail en le 
généralisant et en précisant certains des points qui y 
sont traités. Nous examinerons ici le problème sous 
sa forme la plus générale, réservant pour un article 
ultérieur l’application à la détermination des vies 


moyennes. 


D'une façon générale, nous poserons le problème 
de la façon suivante. Soit x une certaine grandeur 
mesurable pour des particules individuelles ; suppo- 
sons que la probabilité pour qu’une particule corres- 
ponde à x compris entre x, et x, soit une fonction 
connue W (Ty; Los Os Los +. OI) À Li, À et de l para- 


mètres œ, %, ..., œ. On observe un ensemble de 


particules parmi lesquelles n, ont un x compris 
entrer; et'%, 1 entre-m, et... Nr entre YA 
et x, ..., et l’on se propose de déterminer, d’une 
part, la valeur la plus probable des paramètres «;; 
d’autre part, les écarts quadratiques moyens do, 
c'est-à-dire la précision avec laquelle les &; sont 
déterminés. Dans les problèmes traités par Peierls, 
la grandeur x est le temps, mais ce pourra être tout 
autre chose, par exemple l’énergie. 
Nous poserons 


CAEN M ne A Ac 


= DUT, MS Es dan) 


D=— = (ak Œo 84 41) 
V 
Ex 

m = > nie 
V 


Les «, étant supposés connus, la probabilité pour 
qu’on observe précisément n, particules dans l’inter- 


(LS NTo Gilet 


(2) Proc. Roy. Soc., London, Série A, avril 1935, n° 868, 
vol. 149, p. 467-486. 


valle (tot), dans lmtervalle (52): %est 


donnée par 
M! ï 
QE (RC SU UN EN ET 0 
1 ra ne ON EM nn Eee #" ) 


M étant le nombre total de particules qu’on aurait 
observées si toutes les valeurs possibles de x avaient 
été réalisées dans les expériences. 

Admettons maintenant, comme on lé fait cou- 
ramment en probabilité des causes, que la proba- 
bilité a priori pour que M ait une valeur déterminée, 
supérieure à m, est indépendante de cette valeur, 


c’est-à-dire que toutes les valeurs sont également 


probables, et que la probabilité a priori pour que «: 
ait une valeur comprise entre œ; et œ + dx; est 
proportionnelle à da«;. La formule de Bayes nous 
donnera la probabilité a posteriori 


p da dus... day 


ar fau. da D P 


Frs 


Le dénominateur ne dépendra évidemment plus 
de M ni des «;, mais seulement des n, et l’on pourra 
écrire 


SE) M 
= A EN VASE JON TR AT CET OR TA FEI ST ee 
Caub rlnel...nx! (M—m)! 
x mis.  PER(I — JR dos... duy (1) 

ou encore 

Re 17! 

PEKIN (LD): 
\ AUTO ve 


m ! pit sas (ES) 


da... duy. 


AS EN El! Tour 


Si M et m sont tous deux très grands, un calcul 
classique donne 
M! (n— Mi)? 


TA 
2 : ent (1 — CON METT Ye EL EN MeEuT "A EURE M (Is) ; 
m'(M—m)! ) Ten 


RDC AC E FE 


ÿ NC y. ‘ 
Nous pouvons en conclure que P n’a de valeur 
appréciable que si m 72 Mw; or, nous ferons constam- 


ment l’hypothèse que les résultats expérimentaux son. 


ceux qui rendent P maximum. Donc 


TS PUIS 
D = PPS" ee 


K(n,) m | 
w) ulnol...ngl 


Var m(i— 


cm 


En négligeant le facteur Vi—w, qui varie peu, 


on écrira encore 
log P = const. +Y n, logæ,— m log. (2) 


V 


Les paramètres «; devront alors être tels qu’ils 
rendent P maximum, ce qui conduit aux équations 


à ny dw, 
di % 


y. da; 
V 


m dw 
L4 Où; 


4 


en nombre suffisant pour déterminer les œ;. 
- _ Pour une variation assez faible dœ; des «&;, la proba- 
bilité P deviendra 


a log P 
o P — Ôa; da : 
log P = logP ETS 2 are ne ) Be, 


d log? É 
le symbole SR » que nous représenterons 
È 0 * 

" ARRETE ! 9 log P 
aussi par À;; signifiant qu’on a remplacé, dans RE nn 
CIOEET 


Tes œ par leur valeur tirée de (3). La somme double 
qui figure dans le second membre est une forme 
quadratique en dx; qui doit être définie négative 
pour que P, soit effectivement un maximum. Une 
substitution linéaire HO nbE permettra de la 
mettre sous la forme — ZÈ,f$*. Les valeurs quadra- 
tiques moyennes sont 


(EE), 2 


d’où l’on peut déduire les valeurs moyennes dx du. 
On peut effectuer le calcul très simplement de la 
facon suivante. Isolant un indice kX bien déterminé, 
on a 


un Day = Axx dai + 9 Ÿ Aix da Ô pe + 2 tee 


LEk 


( = 
Ps Ge = O0, 


da — 
PS M; 


: 


— VA ze dax + 


+Y A;ij da; da j — 


ijZk 


da, ne figurant que dans le premier terme, on est 
sûr que 
DE 


er 


da; \ dax VAu+ — 
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ou | i 

DA San = 0 (iHR), 

i , 

DT da ; 
ss jZk 
dax VA as = — 1 
VAr 

ou 


(Zn Med 
NA 


Or, l’on peut écrire 


(Zu) . ja day QT ur 


On a donc, d’après les deux relations précédentes 


4 CCR per MALE 
— Ai= Arr A ja 5aj dar ou DA ja dej dot =—7. 
J Î 


F 


Tous ces résultats peuvent être groupés en 


S S Ftenel U Sr k 
DÉTELTEECTES k ne » 

ns 
et ces relations sont vraies pour tous les indices à 
et k. Donnant alors à i une valeur fixe et faisant 
varier X, nous obtiendrons / équations linéaires avec 
l inconnues dx; d&;. La résolution est facile. Soit À 
; o0A 
04; 


le déterminant des A;; et B;; — ; on à 


DHNÔR) A LA (4) 


ce qui résout entièrement le problème. Si par 
exemple il n’y a que deux paramètres « et &, On a 


IN — A1 A5 — A?;, 


On peut remarquer que, dx? et du: étant certai- 


nement positifs, cela implique 
Au < 0, Aa 0, 


et l’on en conclut 


da > — 


Au 


ce qui donne une limite inférieure de dx? et da. 
Le problème important est donc le calcul effectif 


des A;;. On a 
0? log P Ÿ ny 22%, > ny, dw, dw, 
da; da j 5 [Lan d®; da; y, dt; da; 
Y Y 
mm dwæ, m dæ dæ 
æ da;da;  w? du; da; 


\ 


Mais pour calculer À;; — ( er 


da; da; 
Ÿ j 
sans grande erreur, remplacer n, par sa valeur 
moyenne Mw, et m par Mw et se servir de 


ve x < 4 “ 


Il peut se faire que les diverses grandeurs soient 


Ÿ Ales ns, Fée o 
ÿ | ñ L 


it 
\ 


Dan) Van 4 
2 : @; Ôa} 2 ME + D AB. 


Résolues comme précédemment 
. donneront 


TU 


ces équations 


(Sao) = VS auBh 
pou UE AU NE AS 


1 AL 
= Rd AuBxB 
Fr 


on peut, 


DE HN Pwy :  dw 
rot Jay da; — 0a;0aÿ 
x V »/ 
d’où 
ve GENS pr 2, MAT 0? æ 1,8 92 
de dar da; da; da; 7 da; da; ww dada; 
REV Ÿ 
et É | | ; 
“AVE m dœ 0w ny dy dy 
"71? da; da; re da; - 04; 
f { } | sy : 
Ai; I dw dœ nos du, AURA 
.M w\ da; da; D y) da; da; ! fe 


È DA : 
à RU | 
Cas particulier. — Supposons que les inter- | 
‘valles t,—%x,_; soient PAU petits et Cu 
par suite, : 4% 4 


w, = w' (x) dx. 


- 


sf eee [car æ(xo) = 0], 


et prenons comme variable 0, l'extrémité x elle- 


sothssnissessisécisien pénal 


même. On a EE 1 
Air dm dm [1 0#(r)dw (Er) 
M  w dax du . (x) dar day ï 

: t : ! s M 
d’où. , TR 


fonctions d’un paramètre 0 autre que les «; et que En MUR ET OI A C1 dw (zx) 0% 
lon veuille calculer les dérivées de 0 da; par |A dax où POSER idep ed 
rapport à 0; nous désignerons ces dérivées par un 1.0w Qw(æe) 1. dw(zr) dx) à 
accent. Partons de : æ dax dus. CCS ANR TEE AE | 
; 1 dm 1 dw(xx) de * 
DOPTEUT da da; =— Le re É dax  W(xk) Dax ] ; À 
: 1 dw 1: 1 0w'(TE) 5 i 
ee ads | 
Donc | 


EDS, es 


“ 


Comme w" (1) est certainement positif, on voit 


que (da?) est toujours négatif, c ’est-à-dire que da? 
est une fonction décroissante de æ.. La plus grande 
précision s’obtiendra donc en prenant Tr aussi 
grand que possible, contrairement à un résultat 
indiqué, dans un cas particulier, par Peierls à la 
suite d’un calcul numérique un peu approximatif. 


. Manuscrit reçu le 28 juillet 19474 


X 


ner is nnle Suis st hist mat siDos so 


à 


dates d don pl Se à 


1 
U 


THÉORIE DU CÂBLE COAXIAL. I. 


Par F. POLLACZEK. 


Sommaire. — L'auteur établit, au moyen des notions el méthodes classiques de la Physique mathé- 
matique (fonctions de Green, équations intégrales, développements en séries suivant des fonctions 
propres), une théorie rigoureuse du câble coaxial en régime périodique qui est fondée immédia- 
tement sur les équations de Maxwell. Il étudie particulièrement la partie du champ qui ne subit pas 
d’affaiblissement exponentiel et qui échappe complètement à la théorie élémentaire des câbles fondée 


sur les équations de Kirchhof, 


1. Introduction. —— Dans ce qui suit, nous 
esquisserons, pour le câble coaxial droit, de 
section circulaire, une théorie rigoureuse qui est 
. uniquement fondée sur les équations et la théorie 
de Maxwell et non sur les notions secondaires de 
résistance, inductivité et capacité. Nous supposons 
que le câble, de longueur finie ou infinie, est plongé 
dans un milieu homogène infiniment étendu en 
direction radiale et qu’il consiste en deux conduc- 
teurs séparés par un isolant homogène; cependant, 
nos formules s'appliquent dans leur partie essentielle 
à un nombre arbitraire n de couches homogènes 
coaxiales, caractérisées respectivement par les valeurs 
de leur rayon extérieur r;, leur constante diélec- 
trique «;, leur perméabilité pu; et leur conductibi- 
LE EN CNE SAN LE 

On sait que dans bien des domaines des Sciences, 
on est amené à étudier les fonctions propres et les 
“valeurs propres de certaines équations différen- 
tielles ordinaires, ou à dérivées partielles, sous des 
conditions aux limites homogènes; on définit des 
fonctions de Green correspondant à ces problèmes, 
on établit les équations intégrales auxquelles satis- 
font les fonctions propres et grâce à des méthodes 
devenues classiques, l’on parvient, à résoudre les 
divers problèmes aux limites qui se posent dans les 
Applications. 

Tous ces problèmes et notions ont des analogies 
dans la théorie du câble où, toutefois, l’on aboutit 
à des équations différentielles et à des conditions 
aux limites où figurent des fonctions données 
complexes; mais, jusqu’à présent, une seule de ces 
analogies, celle concernant les constantes de propa- 


gation du câble et les champs qui leur correspondent, : 


paraît avoir été étudiée, bien que de manière incom- 
. plète. ns 

: Nous verrons que le spectre des constantes de 
propagation du câble coaxial que nous repré- 
sentons dans le plan de la variable complexe y, 
est constitué par un ensemble (non nécessairement 
infini) de points + Yy, Æ y» ... (spectre ponctuel) 


et en outre, par deux courbes qui sont symétriques 


par rapport à l’origire de ce plan (spectre continu). : 


En des points situés à grande distance r de l’axe 
du câble, les composantes des champs propres de 
première espèce qui correspondent au spectre 


e—const. 7" À 

———, tandis que les 
: 

composantes des champs propres de deuxième 


ponctuel, sont de l'ordre de 


Éore > 
espèce s’y comportent comme “ Deux champs 


propres différents sont orthogonaux en un sens qui 
sera défini plus tard et cette propriété permet de 
résoudre les problèmes aux limites les plus simples 
(problèmes de Dirichlet) au moyen de dévelop- 
pements de Fourier. 

Pour des valeurs moyennes de la coordonnée 
axiale z, mesurée depuis la source d’énergie, le champ 
du câble est approximativement égal au champ 
propre de la constante y, de partie réelle minimum, 
et proportionnel à e-%*, Mais, pour z tendant vers 
l'infini, le champ du câble est égal, à une fraction 


négligeable près, au champ qui correspond au premier 


iOz 


e (5 


z2 


élément du spectre continu, et proportionnel à , 


où c désigne la vitesse de la lumière et « la pulsa- 
tion de la tension appliquée. 

Les méthodes que nous appliquerons ont, pour 
une grande partie été exposées dans un Mémoire 
antérieur (?) qui ne sera pas supposé ici comme 
connu; cependant, dans le présent article où la 
propagation dans des systèmes à symétrie de révo- 
lution est étudiée de manière rigoureuse et plus 
détaillée, une partie de la théorie a été développée 
de manière différente et plus simple que dans le 
Mémoire mentionné. 

Pour être bref, nous ne considérons ici que des 
champs qui sont sinusoïdaux, de pulsation donnée ©, 


(1, F. PczLAczEeKk, La propagation des phénomènes élec- 
tromagnétiques dans des milieux à structure cylindrique 
R. G. E., 1936, 40, p. 103-112 et 131-147. 


A6. 
par rapport au temps et nous supposons, en outre, 


RSS - d 
que ces champs sont de symétrie circulaire (£ = o) 


et que leur composante magnétique axiale H, est 
nulle; d’ailleurs, la théorie des champs que nous 
excluons ainsi est tout analogue à celle des champs 
_ de symétrie circulaire. 

Nous utilisons des unités de Gauss; toutefois, 
vers la fin de cet article, certaines applications de 
nos calculs seront données en des unités électro- 
magnétiques pratiques. Les composantes de tous 
les vecteurs qui figurent ici sont des fonctions 
complexes des coordonnées cylindriques r et z et 
de certains paramètres réels ou complexes, multi- 
pliées [sauf pour EP), H°, E(y), E,, H(7), H,,] par 
le-facteur 2/27, 


2. Spectre des constantes de propagation et 
champs propres. — Pour des états sinusoïdaux 
en {, donc pour 


— = 16, 


il 
ot K ) 
les équations de Maxwell prennent, dans un milie 
homogène, la forme 
ro ARS + RE 


PIC 


top 


rotE = — H, 


C 


où l’on posera 


E — €; eu.) Ex .. TC — O1, +... O4) 


selon les quatre milieux de la figure 1 : conducteur 


, Fig. 1. — Section normale du câble coaxial, 


intérieur, isolant, conducteur extérieur, diélectrique 
extérieur. Nous utiliserons quelquefois la forme 


intégrale de la r'e équation de Maxwell, à savoir 


fa as= FAT (2) 
C C 


où Î désigne le courant qui passe à travers une 
surface bordée par la courbe C. 


En multipliant la 1Te équation par — SE et en 
posant 
ss" À LO + »2 
7 A a (AT po ER OO) 


C2? 
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ON. à 


ADNI LU = 
J rotE — — ÉH. 


rot H — RE, 


On sait qu'avec les hypothèses H; = 0 et. 


De 0, (5) 


les composantes H, et E, s’annulent elles aussi. 
Il vient donc ; 


H.=H,=E£,=0, (6) 
| 


de sorte que les équations vectorielles (4) sont 
équivalentes aux trois équations scalaires 


MAUR ina PM et den 
-r or Je . 
al 
AUS H. 26 FE; m2 dE; 
MAN Ôr 


Nous cherchons d’abord les solutions partielles de 


ces équations qui ne contiennent la variable z que” 


par un facteur e-Y=; ‘en outre, sur les surfaces de 
discontinuité r = r;, r», T2, leS Composantes tangen- 
tielles E; et H, doivent être continues selon les 


postulats de la théorie de Maxwell. Enfin, en appli- 


quant cette théorie au cas d’un câble de longueur 
infinie, nous ne considérerons comme physiquement 
réalisables que des champs dont les composantes 
disparaissent, pour r tendant vers l'infini, au moins 


I 
comme = 
En posant 
E. = et: EN(r), 
E,= elt-ts EQ(r), (8) : 
He = ewt-Ys HU(r), 
nous tirons des équations (7) 
ae Toy ,» tou 7 d£Z) 
BR EN = — AS HUE ee HN= — y EN de 
et avec la notation 
= + y? (Æ HÉNPOUL CE MEC A AN) NE 
nous obtenons 
+ dE | cl? dE 
FAN res HE Le 
À #2 dr ? toux? dr (0 


Pour fixer le signe de x, convenons de désigner 
ici par x, sauf pour kÀ + %?> o, celle des deux 
racines carrées du 2° membre de (9) dont la partie 
imaginaire est positive, donc 


Im(x)> 0. (11) 


De même, nous comprendrons sous X [voir (3)] un 
nombre à parue imaginaire non négative; mais, 
comme alors k tend vers un nombre négatif si la 


conduetibitité tend vers ZÉTO, nous aurons 


SET PANTEU u) va 


pour ç6—0. (12) 
En or. les expressions (8) et (10) dans la 
ire équation (7), nous obtenons l'équation 


. r ke? d£ÆM re EN 
APTE | + n JU = 0 


(13) 


qui, de même que toutes les formules précédentes, 
est valable pour des 6, pm, o qui sont des fonctions 
continues quelconques de r. 

Avec les notations 


EM(r) = uy(r), ie = pr} (14) 
(13) prend la forme 
d u’ 
a (en) (15) 
et cette équation devient 
+ Eu + x2u = 0 (= A+ 72) (16) 


si €, p, os sont constants par morceaux; c’est la 
forme simple que prend l'équation des ondes à 
laquelle chaque composante rectangulaire de E et H 
satisfait dans un milieu homogène, dans notre cas 
particulier où 

D) 0) d 


— —=10 — —=—) 


dt oz ds d9 


=. 0. 


La solution partielle selon (8) doit satisfaire à 
l'équation (16), auf pour r = Lis To To points des 
discontinuité de k et u; mais en ces points, les 
composantes tangentielles EM = u, et HŸ doivent 
rester continues. En outre, toutes les composantes 


« : . . I 
doivent s'évanouir, pour r —> 00, au moins comme = - 


En désignant par Or‘) des grandeurs «(r) 
telles que 
Fr Cre pour 7 
(ou, le cas échéant, pour r +0), 


on a donc [voir (10) et (14)] pour un champ physi- 
. quement réalisable, outre (16), les conditions supplé- 
mentaires suivantes : 


du : 
uw et p— sont continus 


dr (17 


donc pour or <, 


o(:) pour 7->œ. (18) 


pour 


atr)=0 (5); Le 


L’équation (16) ou 


= F1, lo, 73, 


+ +y = (19) 


est l'équation des fonctions cylindriques d’ordre zéro. 
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- Rappelons que (19) possède comme solutions linéai- 


rement indépendantes les fonctions 
d'ordre zéro, de rre et 2° espèces, 


de Bessel 


vly! 
0 


D DE 


MT 


USENET ONE À) 


NZ) — SJ (æ)log EE 


ylyl 


et que, pour les grandes valeurs de |x|, on utilise 


- les fonctions de Hankel d’ordre zéro, de rre et 


29 espèces : 
H4)(æ) = Jo(æ) + iNo(æ), (21) 


HE (æ) = Jo(x) —iNo(x), 
Dans nos formules figureront aussi les dérivées de 


ces fonctions qui sont, au signe près, les fonctions 
cylindriques d’ordre 1, donc 


dJ. 
Ji(æ) = — . 4 
(1) 
HO EE as (22) 
(2) 
Ho (2) = CE), 


Pour les fonctions de Hankel, les formules asympto- 
tiques suivantes sont valables dans le plan des x. 
découpé le long de son demi-axe négatif 


HE (x) = 


suro(s)] 
TiX “A 
VE 
—ix ÿ 93 
HE(æ) = e—ix ji" H+o(i)} (23) 
R'ELe 1e 
= 
QUE 0 ARE) pour |æ|[—. 


Nous utiliserons, en outre, la formule bien connue : 


| Hz) Hÿ(æ) | 4 (24) 
H(x) H®(x) 7x 

Comme la solution u,— El de (16), qui n’est 
définie qu’à un facteur constant près, est finie à 
l’origine, on peut poser u, = J,(«r) pour or ZT, 
tandis que, dans les autres intervalles, u, peut être 
exprimé comme forme linéaire en Hf” et H. 
On a donc, pour certaines valeurs de y qui seront 
précisées dans ce qui suit, 


Jour), OT ST 
«aH$) (rer) + @&HS)(wr), rr£rs, 95 

FO Len as HO Ce r)+ Her), 72ZLrZLr3, (s 
as HO Gur)+ HS) (ur), r3Zr<o, 


, &ç Sont des constantes. A l'infini, u, doit 
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À 


A l'aide des formules : 
à à 
Hvky ALES | 
la condition (17) de continuité des fonctions 4 . 
et Le pour r =; I» Ta S'eXprime au moyen des 


être nul; mais comme selon (23), 


HP Cur)| = æ pour Tr 5% ice 


Im(u)> 0, 


SRE 5 \ ù k 
on a nécessairement, pour la constante «a, qui 
figure dans la 4e équation (25), 


dÿ=0, sauf pour x,20. (26) six équations linéaires et homogènes suivantes : | 

| 

Jo(airs)— @ Hf(wri) — a H@)(rir1) | " AO || 
1014 (ua ra) — di Pa HU (42 ra) — 262 HP) (tra) = 0, | 
| di H$(rrs) + > HE) (rs) — as HW (asr2) — du H° (23 ro) SE (28 

, do H9 (4 75) + do Vo HP) (% Ta) — 3 °3 HO) (Ce To) — dy, Pa HP) (Ce To) KA 0, > | 

as  HW)(rr:) + a HP) Gs 73) — à HS) Gur:)— a; HT) = 0, Ÿ° | 

1:94 

4! 


Œ: Pa H! UC ra) Fa di Va HP) (es T3) m7 147 MC 1) Te di Va HP (C7: r;) 0e 


Supposons d’abord, selon (11), que Im(k,)> 0; alors, la constante a; est nulle et les équations (28) . 


ne peuvent être résolues que si le déterminant 


Jo(xari) HS) (x 71) HS) (ur) 0 0 0 | 
Pa Ja(riri) pe HS (ri) vaHP)(wr1) 0. 0 0 
DXYE AS Hi Cr) HW)Care) HW (xsra)  HS/(xsr2) 0 (49) * 
0 HU(un) MH) vs HW(ur) 0H ere) 0 
0 HO n HO (zsrs)  H@(asrs)  HS(uT3) 
0 o 0) es HY)(xsr3) PH) Crsrs) 0H) (Krs) 


 S’annule. Tous les y susceptibles d’être des constantes 
de propagation de notre câble doivent donc satis- 
faire ou à l’équation 


none — ini er 


D(r)=0; (31) 
nous excluons par une hypothèse les cas où, pour 
certains y, ces deux équations sont satisfaites. simul- 
tanément, cas dont l’étude n’apporterait rien d’essen- 
tiellement nouveau. | 

À l’aide des formules (20), (21), (22), on voit 
qu'aux points y = ik, et y = Hik,; où respec- 
tivement x et x; sont nuls, la fonction analy- 

tique D (y) est régulière: par contre, en y =Æik;, 
D (y) a des points de ramification logarithmique. 
Mais, comme l’on doit avoir nécessairement 


NRA 


_ (30) 


ou à l’équation 


Im(ri) > 0, 


afin que la fonction u, (r) selon (25), (26) s’annule 
à l'infini, nous n’aurons à considérer D (y) que dans 
le domaine où cette inégalité est satisfaite en vertu 
de nos conventions antérieures, donc dans le plan 
des y, découpé le long des deux branches de la 
courbe (30) qui est représentée figure 2. Dans ce 
domaine, D (y) est une fonction paire de y et l’équa- 
tion (31) y possède des solutions + y, (» = 1, ...) 
en nombre fini ou infini qui sont les valeurs propres 
de l’équation (16) dans les conditions (17), (18). 
Car pour ces y,, les équations (28), (26) peuvent 
être résolues pour les a; et les fonctions propres 


correspondantes ne — El(r) selon (25), (26) four- | 
nissent, à l’aide de (8) et (10) des solutions partielles 
des ÉQUARONS: du câble. 


Plan des Y 


Fig. 21=—"La "courbe Im (x) 0: 


Comme nous l’avons indiqué dans l'introduction, ! 
ordinairement seules les racines +%, de partie |! 
réelle minimum sont prises en considération pour 
la propagation de signaux téléphoniques ou télé- | 
graphiques. Pour les “valeurs usuelles de 1, rs, rfi 
et de w, les nombres rixf) = r, VK? + y? et rx), à || 
i—1, 2, 3 sont petits, si bien que, dans (28) | 
et (29), les fonctions H{° (x), ..., H° (urs) se 
réduisent à des expressions simples; en outre, en 


ANT TH 
raison de là petitesse de y, 4 et x) peuvent 
respectivement être remplacés par k, et k,. Dans ces 
-Circonstances, l'équation (31) fournit, pour +, la 
même valeur, et les équations (25) à (28) donnent 
pour le champ les mêmes expressions que la théorie 
usuelle des câbles. | 
Outre les solutions propres ordinaires ou de 
17e espèce u, (r) dont les valeurs propres + ET 
forment le spectre ponctuel de (16) [avec (17), (18)], 
il existe une variété continue de solutions propres 
de 2e espèce u, (r) correspondant aux points y de 
la courbe (30) (fig. 2) qui forment le spectre continu 
de (16) [avec (17)]. On obtiendra uy(r) en résol- 
vant d’abord les deux 1res équations (28) pour 4, 
et a, ensuite les 3e et 4e équations pour 4, et a,, etc.: 
car les déterminants [voir (24), (27)] 
4 ik? 


d, D RSR - = 
Ty l'y 12) 


(y as, 4) (32) 
de ces paires d'équations sont différents de zéro. 


En particulier, on déduit ainsi les formules 


a D*(x) Les D(Y) 90 
d; d: de, HE dd; d,, (ES. 
où 
k? r 
D (x) = rt Gurs) A6 (1) + HG rs) Ase(Ÿ), 
, (34) 


D*(1)=— SR Her) Ass (7H Xasrs) A50(Y), 
A 5% et A6 désignant ici des mineurs de la 6e colonne 
de la matrice (29); la 2° équation (33) met en évidence 
que pour les points y, du spectre ponctuel, u, se 
confond avec u,. 

On voit aussi que la fonction u, (r) formée ainsi 
est une fonction entière de y, ce que nous utili- 
serons par la suite; mais il faut retenir que u, (r) 
n’a pas de signification physique, sauf pour les 
points des spectres ponctuel (31) et continu (30). 

Il nous reste encore à étudier les fonctions propres 
de 2° espèce pour r oo. Or, les formules (23) 
montrent que pour x, réel et r — co, le 2° membre 


de la 4° équation (25) est seulement of): 
; & 


les u, (r) de 2° espèce ne remplissent donc pas la 
condition (18). 

Mais en intégrant e-Y*u,(z) sur un arc quelconque 
du spectre continu, on obtient une fonction E:(r, 2) 
dont on tire une solution E:, E,, H, des équa- 
tions (7) du câble qui satisfait aussi à toutes les 
conditions supplémentaires de la théorie de Maxwell. 
Pour montrer cela, formons, avec une fonction 
analytique f(y), les intégrales 


V2 
Er, 2) = ‘ eiot-yzu(r) f() dr, 

Y1 

Ve ; y l À 
(na =—f eme Eu(rf(nan à (85) 


it 


Ye : cl? ; 
Je etais Es &(r) f(x) dv 
Y1 > 
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duy(r) 


où u,(r) =. r—; ici, la sommation sera étendue à 
ar 


CAD 
un arc quelconque y,Y, de la courbe (30), dont les 
extrémités sont supposées différentes de + ik,. Ces 
fonctions satisfont évidemment aux équations (7) et 
aux conditions de continuité requises pour r —r;, la; las 
car u,(r)satisfait à (16) et (17). Pour démontrer qu’en 
outre, les fonctions (35) sont au moins de l’ordre 


ll . 
de . pour r->00, nous y remplaçons u; par le 


2e membre de la 4e équation (25), et y par la variable 
d'intégration réelle x, —x; on a alors des inté- 
grales de la forme 


A 
lue [| H(xr)e(x)de  (u=0,1 v=1, 2), + 
x" \ 
où 4 (x) désigne une fonction complexe de x. 
Or, à l’aide de (23) et en utilisant le fait qu'ici, 
beton 4par exemple “pour L'— 0 VIT. 


LA ti D 
= |. ie Rs of Es 
rat y VX rVr 


Loi = 


2 ‘ 


La Ê * gx) irx ts 
= cie - d(e ose) 


; Ty 
tr — 
2 


ae 
S 
ele 
à 
Æ 
=. 
me, 
ES 
-G 
RS 
SIIR 
> 
@ 
= 
[RE | 
> 


- Donc, les fonctions (35) forment une solution 

légitime des équations (7); en y remplaçant f (y) 
AUOT DU : 

PAT 55 Dr) [voir (33)] { en aan tendre le 

point y, vers —ik,, on obtiendrait même des solu- 

tions qui sont exactement de l’ordre de . pour r —>00 


et les développements ultérieurs montrent qu’inver- 
sement, on peut représenter chaque champ maxwel- 
lien à l’aide des champs propres des deux espèces. 


3. Orthogonalité des fonctions propres. — 


Nous démontrerons maintenant que les  déri- 
n d * ï 

vées roi et de deux solutions propres 
7e ; 


différentes sont orthogonales. Pour cette démons- 
tration, nous nous servirons de l’équation (15) qui 
est plus générale que (16) et n'’utiliserons pas 
que k? et y sont constants par morceaux; par consé- 
quent, les formules qui suivent seront valables, 
du moins formellement, pour des k(r) et des m(r) 


qui sont continus par morceaux, et constants à 
partir d’un certain r = R;. 
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Soient maintenant uw, (r) et &,(r) deux solutions 
propres de (15), appartenant aux valeurs , 
et ÿ7<+ 7, du spectre ponctuel, donc des fonctions 
qui satisfont, la première pour y — y, et la deuxième 
pour y = y, aux conditions (15), (17), (18); on a alors 
la formule 


u( (Gal) mie 
ref p( Dir D Re dr = 0 


[r- É EYi p(r)= |: 


En effet, on obtient à l’aide de deux intégrations 
partielles qui sont légitimes en raison de (17), 
Wir) ue (Tr) “+ 


R 
Di PU) LRU EYE 


(37) 


Îl 


d PU ae 
. LL )| dr (38) 


et comme la fonction à intégrer du dernier terme 
est nulle en vertu de (15), il vient 


a pUuiUs  pusus JR 

ICR) = = |rerr pe | 

ULRR)  [ui(Rhu(R) ME (39) 
Mn ra US So RENE 


mais grâce à (18), le dernier momie tend vers zéro 

pour À — co, ce qui démontre la formule (37). 
Remplaçons dans ces équations w,(r) par une 

fonction propre de 2° espèce u, (r); dans nos hypo- 


a) HO (rx) + XL) HUE) 


2 Te 
und. pi | Cat) AD G) 2 007) HE) Le (D AUD (ren) + a) AN | TE 


Il faut maintenant distinguer deux cas. Admet- 
tons d’abord que ni le point y ni — 7 ne sont situés 


JUT 


ACT 
sur l’arc y, y», de sorte que le facteur - sy est 


borné; en appliquant alors la formule (36), on voit 
que (44), donc aussi le re° membre de (43), est nul. 


Supposons maintenant que le point y est situé 


ANG 
sur y,7y, et traçons dans un voisinage de y deux 
arcs C' et C” qui soient convexes respectivement 
vers le haut et vers le bas (fig. 3). Admettons que x; 
est positif [ce qui est le cas, en raison de nos conven- 
tions sur le signe de la racine de (9), respectivement 
sur les bords supérieur et inférieur des branches 
droite et gauche de la courbe (30)]; on a alors, 


EURE Im(x,) Lo, 


thèses, l'équation (15) se confond pour = :| 
avec (16), de sorte que l’on a [voir la 4e équaz N 
tion (25) et les formules (23)] : 


du Cr) = a) Her) + 0 (DH Gr) (0) 
pour, Im(xi)=0, re À 


uy(r) = 0 () u(r) = dt (M0 
F \r ra 


pour 7'—> cc. 


Dans ces conditions, Z(R) (38) tend encore … 
vers zéro pour À -> 00, de sorte que la formule (37) 
reste valable si l’une des deux fonctions propres | 
estidera° espèce. 

Mais si u, et u, sont tous deux remplacés par des 
fonctions propres de 2° espèce, la condition (4r). 
ne suffit pas pour faire disparaître I (R) pour R — co. 
Cependant, dans ce cas, la formule (37) subsiste 
entre les fonctions u, (r) — u, (r) et ù 


u2(r) = — = 2 . k: uy(r) FC) dy (42) 


Ÿ1 


f(») désignant une fonction analytique et l’inté- 
gration étant étendue, comme dans (35) et (36), 


- AIDES 0 
à un arc y,, y, du spectre continu. En effectuant 
dans l’équation (9), AVEC W =U,, U=Uy, l’opé- | 


ration lim  . .f(Y) dy’, on obtient 


Le P(r) D i . f(x") dy’ us — dr 
DA À 


2x). ‘K?+r 2+ y 
2 
Y: uy(r) un ft) IA 
I DR ; | 
= lim ea pr) | uy(r) u,(r) JO qy (43) |] 
ro 271 A DEL a Y?— y! 


1 ï 


et à l’aide de (40), on a, pour le 2° membre, l’expres- | | 
sion | 


at (y) Ar) + ay) HET x) 


à l’aide de (23), 


sur C: Im) > 0, Et 
e—/ 1m (x) 


= pour 7 —c,.1|| 
Vr | 


| HD (4, r) PA | HE r) | vo 


r'Im(xi) 3 
[HDi r) lo LH (x) | © He pour 700.4 


‘| 
té ñ 


En remplaçant le chemin d’intégration de (44) 
dans un voisinage de y par C' ce qui ne change pas | 
la valeur de cette intégrale, on obtient done que | 
les termes qui sont multipliés par HU(x;r) ou .| 
par H (x; r), donnent pour résultat zéro; en rem- 


1 
plaçant ensuite, pour les autrestermes, St A dy’ 


a 7 
par — sy oi CSC SAN) NOoUsS voyons 
que tot  . ÉSal à à la limite de ce résidu, donc à 


[HP Gur). HKGur) f(X) 


lim k jU) )at2) 
re Hi% (M) tr) HP CET A (er) | 2% 
2ouAt AE 
— RES CS (2) 
a CCR 


NE 3. — Déformation du chemin d'intégration de (44). 


On obtient ainsi entre la fonction u, (r) selon (15), 
(17), (40) et la fonction (42) la formule d’ortho- 
gonalité 


a ; Y: u (Cr) LS UC) d Rp 
T) — Y )dy : r 2) 
TA RS RSR CO Ar tranee (45) 
û fa l i 
! s TS 
0 SRE EVA, 
En CES CR RU LE D 4 TS 
— + qu) (y ar} 00) si + E vi Y3. 
> , NE Dr i TA 2 
Ti ik d 


4. Introduction d'une fonction de Green G*.. 


— Nous étudions ensuite la fonction de Green de 


" FUN, ve 
Ve RL à Ë ÿ 
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l'équation (16) sous les conditions (17) et (18). 
C’est une fonction G (r, r') de la variable réelle r 
et des paramètres r' (réel) et + (complexes), définie 


pour toutes les valeurs >= 0 de r et r’ et pour tous 
P 


les y qui n’appartiennent pas aux spectres ponctuel 
ou continu; GT satisfait à (16), (17), (18) sauf 

dG) 
au point arbitraire r —r' où nous imposons à QT. 


la discontinuité suivante : 


NENAIRE 2G\W) Ho De ne 
Ê Cris Ÿ? dr EU = {, ( 16) 


Pour calculer G, on procédera pareillement que 


lors du calcul de u,(r) — Et)(r) dans l’hypo- 
thèse (11) qui avait comme suite la relation (26). 
Admettons par exemple que r' appartienne à l’inter- 
valle 

NÉE (47) 


Dans ce cas, G" satisfait à l'équation (16), sauf 
pour r = 73 l', Ta, T3 Ae sorte que l’on obtient les 
formules suivantes : 


(O2 DE 
TANT ON 
bd: Her) + b, HS) (kr), 
ds Hs r) + b, HOT), 
b HS) (ir), 


bo T5(rir), 
bi H) (x r) ag b2H9)(% 7), 
NEA) 


TT Sp 
1aLT< 0, 


dont les coefficients b (r'), dépendent des 
paramètres r' et y; ensuite, les conditions (17) 
ét (46) entraînent [en utilisant les notations (22) 
et (27) huit équations HRÉQATES non homogènes, 
à savoir : 


Lo (ar4) — O1 H\) (22 274) — 0» H? (2274) 0, 
Ja(rars) EE 1) — DH )(rors) nn 
EN HO Gor Nb, HoGer) 6 HU) HE) DEC 

bis H (er) + b50 107") — b:02H (ar) — Dies HET") = 

b: HE Cor») + b HP or) — b; Hasr2) — 06 HS)(rsr2) 108 

D302 H1 toto) + br v21 (asr2)— b:e:H\ sr) me be3H@ 17) = 0, 

D; Hi) Cesrs) + 05 HP (Cesrs) —0 7 Hi Car T3)= 0, 


Le déterminant de ce système est égal au 

produit des déterminants (24) (avec æ —x7) 
et —wD(y) [(29)] et de manière générale, égal 
: tk? 
ä Firxuxir a DO), 
de la figure 1. Les équations (49) et celles qui leur 
correspondent dans le cas où r’ appartient à un 
intervalle différant de (47), possèdent donc une 
solution et G% (r, r’) existe, par conséquent, pourvu 
que D(y) “0, c’est-à-dire que y n'appartienne pas 
au spectre ponctuel de (16). | 

Entre la fonction G% et les fonctions propres 
existent des relations qui sont importantes pour la 


si r' appartient au v‘"° milieu 


b503H4 VCksr: ») + 05031 107 l'a 3) — 0; Cuir D — 0, 


solution des problèmes aux limites et pour le déve- 
loppement d’un champ arbitraire suivant les champs 
propres qui correspondent aux fonctions propres des 
deux espèces. En procédant de manière analogue 
que lors de la démonstration de la formule (37) et 
en tenant compte de (46), on obtient les formules 


ve )= (y? p(r) dGW\(r,r) u(r) dy 
VMS : Je = y2 Dr F4 1 
APE ÿ )) (50) 
à AS 3 cs ARR 
us = ne 72 DRE ET Dr Rey 
(GIE F y), ] 


PT SA AS OC) 


où u, D = — E0) (r) & uy(r) = E) (r) sont respec- 
tivement des de propres de re et de 2e espèces 
de (15) ou (16). 
En appliquant le même procédé: aux fonctions 
G (r, r') et G) (r, r"), on obtient la formule bien 
connue | 
ON OU) (nr) (51) 
On peut se rendre compte facilement que, 
pour y — y,, tous les mineurs A;, de la matrice (29) 
ne peuvent pas être nuls à la fois, ce qui entraîne 
 l’unicité des fonctions propres u, (r) et, de là, on 
conclut que les zéros y, de D(y) sont simples. Dans 
ces conditions, u, (r) est égal, à un facteur près qui 
dépend de r', à la fonction 


lim (y — y») GM(r, r)= 5 f GY\(r, r')dy, (82) 
Y>T PR 


_K, désignant ici un petit cercle autour du point y,; 
_ car la dérivée, par rapport à r, de cette fonction est 
continue pour r —r', de sorte qu’elle satisfait aux 
équations (16), (17), (18), pour y =, En raison 
de (51), on a donc 


lim (7 — 9) Gr, 7°) = œuy(r)uy(r'), (53) 
one 


où c, est une constante que nous pouvons calculer 
en introduisant ici pour r et r' des valeurs particu- 
lières quelconques. 

En choisissant par exemple r — 0, r! —r}, ce qui 
est compatible avec l'hypothèse Un. dans laquelle 
nous avons établi les formules (48), (49) de G, 
on tire de (25) et de la r'e équation (28), 


BP= VRAI 


Ensuite, la r'e équation (48) et les FanatiQRe (49) 
donnent 


uy(o)=1, ar) =Jo[x/ri] 


er 
rD(Y)’ 


où A: (y) désigne le mineur de l'élément dn de la 
matrice (29), et il vient 


GW(o, r1) = 671) = 


An CNE AQU) 
raD(Y)7 raD'(1v) 


En introduisant ces expressions dans (53), on 
obtient 


Him (y —y,) Go, r1) = lim (y —y,) —— 
> TE T 


EU Le Ax(Yv) tv), R 
97 Ada D)? ‘4 
de sorte que (53) prend la forme 
lim (y— y») G(r, r°) LH (Tr) Cr) Ari(yv), (55) 


ST raJo[xP/r1] DC) 


Cette équation montre que la fonction analy- 
tique de y 
uy(r)uy(r) 


r1J6(%171) 


11(Y) 


GT) * / Mi Le 
era D(1) 


(56) 
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est régulière aux points y — y, du spectre ponctuel, | 


quelles que soient les valeurs de r et r'; d’ailleurs, 
on peut démontrer aisément que le numérateur de 
cette expression est algébriquement divisible par D (y). 

Désignons par f* (y) la fonction en laquelle une 


fonction analytique f (y) [qui soit holomorphe dans ! 
le plau y découpé selon la figure 2] se transforme 
quand y contourne le point — ik, une fois dans le 


sens positif; on a, en particulier, 
LH GT = — HP Car) | 
a) 


LUE T2) ASE LE À 
27 Le 297 


et ces formules suffisent pour toutes nos applica- 4 


tions. Considérons maintenant la différence 
AGIT, r)= Gr, r)— G 
entre les valeurs que G' prend sur les bords supé- 
rieur et inférieur de la branche droite de la cou- 


pure (30). Considéré comme fonction de r, AG) 
satisfait évidemment à (16) et (17), même pour r — ie 


et est, d’ après la 4° équation (48), de l’ ordre de — F 
Fr 


pour r—>co; c’est donc une fonction propre de 
2e espèce. Comme c’est, en outre [voir (51)], une 
fonction symétrique en r et r', on a 


AG (7, r qi c(y)u,(r)uy(r) [{m(x) = 0] 


GE reiTE) tree) (58) 


(59) 


et c(y) peut être calculé, comme auparavant c,, 
en introduisant ici pour r et r' des valeurs parti- 
culières. 


Avec les valeurs précédentes r =0, r! EE 
le calcul fait plus haut fournit la relation 
EE A 
raJo(iri) L'D(T) DPCY) 
ie L \ A51(T) 
= ANG) D(1) Se 
et de là, 
Ag UrC) uy (7) À Asi(N), 6 
An) D) ee 
Comme la fonction “(72 07) est pas ramifiée 
rido(x ri) 
en y — +ik,, la dernière équation peut être écrite 
sous la forme 
a fem MCD At Là (ce 
[ere ren) AtD] Lo (62) 
donc, la fonction (56) est régulière aussi bien 


en y — +ik, qu'aux points du spectre ponctuel, 
et cette propriété qui subsiste évidemment pour les 
dérivées (par rapport à r ou à r’) de (56), sera utilisée 
dans la suite lors du développement de la compo- 
sante ÆE;(r,z) d’un champ E quelconque d’après 
les fonctions propres ET) (r) et EM (r). 

Les grandeurs €, (54) et c(7) (60) peuvent être 


réprésentées par des intégrales définies où figurent - 


(57) 


Et RE qu Pains. # # es 


adidas shirt pin sn 7 


les fonctions propres. Pour obtenir la première de 


ces formules, faisons tendre y vers +, dans la 


ire équation (50) en utilisant la formule (53); il 
vient 


à ; a re Le (OT) ii hi 
u UN Cr) 2)CS — \—— | dr u,(r 
v(77) nf æ ES y(7) 


Ÿ 


1, ° A Fu(r)1 È 
& ==; [ ner 1 dr. ; (63) 


A0 \ 


ou 


: 


Appliquons ensuite la formule (45), pour y € un 

à la fonction propre de 2° espèce u, (r) selon (25) 

et posons-y f(y) — y: il faut alors introduire, 
dans (45), pour at = a; et a = a les expres- 
sions (33) et l'on obtient, à l’aide de la notation (32), 


LES tu 
A2 2f — Le = y (r) dy’ LA + 
A 2Té ux? 


T3 D 1" FF) GAS lo 9 lo gd: 
me nr PAT REA G=Æ+ 72). (64) 


On vérifie dont à l’aide des équations C4 
(32), (34) l'identité 
I LATE HD Gxirs) 
DEP. d'A D) | 


A 


grâce à laquelle on a 


æ LC dl CRT) 7. 
À 1 \ = he : AN - Ô . 65 
Lf dr | Ta A66 D(r)> °e 


en outre, on peut établir, au moyen d’un calcul 
simple qui porte sur les mineurs de la matrice (29) 
et que nous laisserons au lecteur, l’identité 


RU Sr) 1 Adult): 
FA DO) | Gain) D() 


(66) 


Des deux dernières équations et de (60), on tire, 
HÉUr (63), la formule 


ie ef Lu (7 ap Ur) AT Vi 
ET prë (67) 


eh (Ge=e+r) 


5. Interprétation physique de GT et des 
formules d’orthogonalité. — Les formules du 
dernier chapitre admettent une interprétation phy- 
sique simple. Identifions d’abord la fonction __ G®) 
avec la fonction EX qui figure dans (8) et (10) 
et formons, selon ces équations, les expressions 


F | 210 \ 
ewiy2 Dr, rl) = eivi-vs aie GNU TS), 
c? 


= 200Y 2GN (7; r' ) 1Q 
eiwt Ÿz E® (7, 5 64 te eimt—Yys en : , : (68) 
e 2 k2 0G fr, 
ott2 HO (r, r) = eivtts 2 PAR AU 
LE aim 
LOT He Cr 7 ) e ee Fe 


\ 


A 
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Ces fonctions sont les composantes du champ 
électromagnétique d’un câble selon la figure 1, 
car G® satisfait aux conditions (16), (17), +» 


pour PT: CU (18). 


Sur la surface latérale de cylindre 
Fer (69) 


la composante H, de ce champ subit, en vertu 
de (46), la discontinuité 


A Hi Leiwt-v: HO (7, Fa ee re _ 
r 


et à celle-ci correspond, d’après l'équation (2), 
un courant de surface 


1 = 2x1 A = eivt-vs. VOL) 

Nous pouvons donc dire que le champ électro- 
magnétique « de Green » (68) est engendré par 
une nappe de courant Sinusoïdal 7 = e”-Y: qui 
s'écoule en direction axiale le long de la surface (69). 
Multiplions maintenant le 2° membre de la 
formule (37) [où u,(r) (14) est l’abréviation 


pour ÆŸ (9) = Er) (n] par 


C5 : C 
2 etDET13( — eiWi—Y:23 9x: 
ÂT ( À Y1) tu) Vy 


compte tenu des équations (8) et (10), (37) prend” 


alors la forme 


= eivt-Vs pl) (r)eivt-V:35 H2)(r)2x7r dr = 0 
jet) Cr) 


(VAE Yi) 


Cette équation exprime le fait que le « flux 
complexe » du vecteur de Poynting mutuel 


A Leon E(r), iv HE2)(r)] (73) 
A7 

de deux champs propres (ordinaires) diflérents à 
travers un plan quelconque perpendiculaire à l’axe, 
est nul; notons que ce flux doit être formé au moyen 
des composantes complexes des deux vecteurs qui 
figurent dans (73). 

Comme des intégrales pareilles à (72) se trouvent 
dans une partie des formules ultérieures, nous 
introduisons, pour le flux du vecteur, 


[EC 2), H(r, r)] 
EA 


à travers le plan z — const. la notation 


ne) ee z)Holr, r')arr dr = E(z).H(r') (74) 
1 


où la variable d'intégration est supprimée; de même, 
nous os 


_- "EU AD (r) arr dr=æElW.H0. (15) 


RERO 


12: (70) 


MR + — 
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Supprimons dans (72) les facteurs eiviys 
‘avec la notation (5), _cétte équation pjens alors 


la forme 


E. H® —EC.H0)—0 Henos (76) 
et, de la Aster on à 

EU), H0) = EU. H0) = 0 (T1) 
pour le fus complexe mutuel de deux champs 
propres de différentes espèces. \ 


Grâce à (74), l'équation (45), te pour des 
solutions propres de (16) qui sont pour r> R, de 
la forme UE devient 


“EU f(r) ay. Hi VUE) 


0 si +yef, 
RU C Me: à ! 
— aft(n) at AAC si 


27Ob 4? rep; 


\ 


c’est'la loi Porthobenalite entre un champ de 


| 2e espèce et rhiee par rapport à dy, d’un tel 
champ. 
Avec les notations (75) et 
1 Ye y nil 
Et)(r) dy’ 
| (79) 
ARON ÿ 
[+ EYrY2; Ha = o], 


les équations (54), (63) et (60), (67) donnent respec- 


X 
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 tivement 


“et, de manière analogue, on tire de . 2° équa- 


Jump. GauUraiER-ViLLARs, 55, quai des Grands-Augustins, Paris. — 1947. ‘ 


ri Jo[ x ra ] D'(ys) | 


Res ln, ao. 
Cy An(Yv) ARR TU M En | 
a Fe Oo 
I 4 mi Ho, Re 
LORD ASE = y EN.H" 
En) 2 D(Y) ce if 


! d'après (68). Pa Eÿ x, : 
les te (53) “ 59) deviennent 


a r) EU). 
L [EG ryar == DGA EE 
271 Jy, ; < t 


SEUL HV 
r') 2 


DÉC) ENNCRES 
11 El). HN) 


AE (7, 


vement par ——— = - et par st? on obtient 


E one ni. | 
VV AE 


EH (7) = 


Yv \ A 
2 — El) | 4 h 2 
dre EU 


+ 


tion (50) es 
Eni").3 He EC | 
’ f ÿ* 4 
El). (Y) = Se ( 
HM (7°) mL EEE Nr). 


(A suivre) 


